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I.  MANAGEMENT  REPORT 


A.  Summary 

We  have  continued  to  improve  our  ZnO  deposition  process  to  the  point 
where  the  films  we  are  producing  are  far  better  than  those  of  which  we  are 
aware.  The  present  sputtering  rate  is  18  urn  per  hour,  a  factor  of  18 
better  than  we  were  obtaining  two  years  ago.  By  carrying  out  deposition  at 
500°  C  ,  we  have  obtained  x-ray  results  which  show  a  full  width  of  half 
maximum  of  the  <0002>  peak  of  0.2°  and  show  very  clearly  that  the  Cu-ka^ 
and  Cu-ko,  lines  are  split.  RED  results  are  of  similar  quality,  and  the 
x-ray  and  electron  diffraction  results  are  by  far  the  best  we  have  every  seen 
and  indicate  that  the  films  must  be  very  nearly  single  crystals.  We  have 
carried  out  further  developments  of  our  Sezawa  mode  convolvers  and 
correlators.  We  have  made  a  convolver  with  an  excellent  efficiency  of 
-50  dBm  with  a  bandwidth  of  25  MHz  .  We  have  made  low  loss  delay  lines 
with  32  MHz  bandwidth  and  have  designed  chirp  transducer  devices  which  are 
presently  being  constructed  with  as  much  as  50  MHz  bandwidth. 

We  have  tested  Sezawa  wave  convolvers  and  correlators  to  determine  the 
effect  of  dispersion.  This  is  in  agreement  with  the  theory  we  have  worked 
out.  As  far  as  correlators  are  concerned,  the  theory  indicates  that 
dispersion  is  far  less  important  than  in  a  convolver,  and  dispersion  should 
not  limit  its  performance  even  in  very  wide  bandwidth  devices.  We  have  found 
certain  problems  in  Sezawa  mode  correlators,  which  we  had  not  anticipated,  but 
for  which  we  have  solutions.  One  is  associated  with  reflection  from  the  top 
plate  due  to  the  much  larger  value  of  Av/v  in  these  devices  than  that  which 
we  had  used  before.  This  problem  is  easily  remedied  by  redesigning  the  top 


plate.  Another  was  associated  with  stronger  excitation  of  bulk  waves  than 
with  the  Rayleigh  wave  devices.  Again,  by  bonding  the  semiconductor  with 
indium  to  a  brass  box,  we  have  eliminated  the  bulk  wave  excitation.  Phase 
distortion  due  to  the  top  plate  has  also  been  a  problem  due  to  the  larger 
value  of  Av/v  ;  again,  with  redesign  of  the  top  plate,  we  can  eliminate  this 
problem. 

We  have  carried  out  a  complete  theory  of  the  storage  correlator  and  can 
predict  all  its  characteristics  as  regards  to  storage  amplitude  as  a  function 
of  input  times  and  amplitudes.  The  theory  is  in  excellent  agreement  with 
experiment;  we  have  written  a  paper  on  the  subject  which  has  been  accepted  for 
publication.  We  have  constructed  the  first  monolithic  Schottky  diode  storage 
correlator;  with  it,  we  have  demonstrated  that  we  can  store  signals  in  as 
short  a  time  as  5  ns  for  as  long  as  1  ms  .  We  have  managed  to  make  far 
better  quality  Schottky  diodes  than  heretofore  by  carrying  out  the  depositions 
in  a  cleaner  system  to  which  we  have  had  access  at  Varian  Associates. 

We  will  continue  with  our  work  on  deposition  of  Si  on  LiNbOj  and  have 
encountered  problems  with  cracking.  We  have  been  examining  a  new  concept,  the 
FET  storage  correlator,  and  are  working  out  a  detailed  design  for  this  device. 

B.  Research  Program  Plan 

We  will  be  constructing  very  broadband  Sezawa  mode  correlators  and 
working  to  improve  our  Schottky  diode  correlators.  It  is  our  intention  after 
constructing  broadband  Sezawa  mode  correlators  to  demonstrate  their  use  in  NDE 
systems. 

C.  Major  Accomplishments 
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We  have  demonstrated  superb  ZnO  thin  film  layers.  The  films  now  look 
virtually  like  single  crystals.  Deposition  rate  is  18  urn  per  hour.  We  have 
constructed  the  first  Sezawa  mode  correlator  and  the  first  monolithic  Schottky 


diode  correlator. 

D.  Problems  Encountered 

No  major  problems  have  been  encountered. 
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II.  TECHNICAL  PROGRESS  REPORTS 


A.  ZnO  Technology  Development 

In  the  last  six  months,  we  have  been  carrying  out  a  series  of  runs  to 
find  the  optimum  conditions  for  sputtering  ZnO  .  Our  efforts  have  been  very 
successful  as  we  presently  have  a  set  of  deposition  parameters  where  the  ZnO 
quality  is  the  best  we  have  done  so  far,  and  probably  better  than  anything 
reported  in  the  open  literature.  The  new  set  of  parameters  are: 
r-f  power  =  1.25  kw 

oxygen  partial  pressure  =  7  urn 

target  to  substrate  spacing  =  4.2  cm 

substrate  temperature  =  500°  C 

sputtering  rate  *  18  ym/hour 


The  films  were  evaluated  by  x-ray  diffraction,  x-ray  rocking  curves,  scanning 
electron  microscopy,  and  reflection  electron  diffraction  (RED).  A  summary  of 
the  results  and  curves  is  in  a  paper  that  was  submitted  for  publication  to 
Applied  Physics  Letters  in  July.  The  x-ray  results  show  a  full  width  at  half 
maximum  of  the  (0002)  peak  of  0.2°  and  show  very  clearly  that  the 
Cu-kctj  and  Cu-ka,,  lines  are  split.  The  RED  results  show  a  very  well 
defined  pattern  of  spots  corresponding  to  diffraction  planes  in  the  film;  this 
indicates  an  excellent  orientation  in  the  films.  The  x-ray  and  electron 
diffraction  results  are  by  far  the  best  we  have  ever  seen  and  indicate  that 
the  films  are  almost  single  crystalline. 

We  will  be  continuing  our  characterization  runs  to  investigate  further 
the  effects  of  substrate  temperature  and  target  to  substrate  spacing.  The 
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ZnO  films  will  still  be  evaluated  in  the  same  fashion,  but  we  will  be  adding 
a  few  more  evaluation  tests,  namely  adhesion  by  indenting  the  films  with  a 
Vickers  indentor,  thickness  uniformity  and  surface  smoothness  with  an  Alpha- 
step  height  profiler,  optical  transmission  of  a  laser  beam,  and  resistivity. 

Measurement  on  resistivity  indicate  outstanding  resistivities  as  high  as 
10^  ohm-cm  .  Preliminary  adhesion  measurements  also  indicate  very  uniform 
adhesion. 

R.  Sezawa  Have  Devices 

In  the  last  progress  report,  we  reported  on  the  fabrication  of  broadband 
(31  MHz)  monolithic  Sezawa  wave  delay  lines.  During  the  last  six  months,  we 
have  fabricated  delay  lines  with  the  same  bandwidth,  but  with  4  dB  lower 
loss  (18  dB)  .  We  have  also  fabricated  high  efficiency  (Fy  =  -50  dBm) 
broadband  (25  MHz)  Sezawa  wave  convolvers  and  demonstrated  them  in  pulse 
compression  experiments.  Finally,  we  have  fabricated  Sezawa  wave  storage 
correlators  and  isolated  several  problems  unique  to  Sezawa  wave  devices. 
Solutions  to  these  problems  are  presented. 

1.  Sezawa  Wave  Delay  Lines 

We  have  fabricated  short  (4  mm)  Sezawa  wave  delay  lines  with  very  low 
loss  (10  dB)  .  Of  this  loss,  6  dB  is  due  to  the  bidirectional ity  of  each 
transducer,  1.5  dB  is  due  to  ZnO  loss*,  .2  dB  is  due  to  resistive  losses 
in  the  1  mil  diameter  bonding  wire  and  tuning  inductor,  and  the  remaining 
2.3  dB  is  probably  due  to  bulk  wave  excitation.  With  the  use  of  the  five 
element  tuning  circuit  shown  in  Fig.  1,  a  3  dB  bandwidth  of  31.5  MHz  was 
achieved  with  18  dB  insertion  loss  (see  Fig.  1). 


2.  Sezawa  Wave  Convolvers 

Sezawa  wave  convolvers  with  an  excellent  efficiency  of  -50  dBm  and  a 
bandwidth  of  25  MHz  (Fig.  2)  were  fabricated.  Only  three  elements  (Fig.  2) 
were  used  for  the  interdigital  transducer  (IDT)  tuning.  The  device  parameters 
are  given  in  Table  I.  The  results  of  5  bit  and  11  bit  Barker  Code 
correlation  are  shown  in  Fig.  3a  and  3b,  respectively. 

A  series  of  chirp  correlation  experiments  were  performed  to  see  if 
dispersion  would  limit  the  usable  time-bandwidth  product  of  these  devices. 

The  experiments  used  3  us  long  chirp  signals  of  variable  bandwidths  up  to  a 
maximum  of  23  MHz  bandwidth.  (A  variation  of  23  MHz  in  3  us  is  the 
maximum  frequency  slewing  rate  obtainable  from  our  chirp  generators.)  To 
avoid  pulse  broadening  due  to  the  nonlinearity  of  the  chirp  frequency 
variation  with  time,  one  chirp  generator  was  used,  and  its  output  was  mixed 
with  two  cw  signals  to  generate  the  necessary  up  and  down  chirp  signals 
(Fig.  4).  Using  the  input  23  MHz  chirp  shown  in  Fig.  5a,  the  compressed 
output  pulse  is  shown  in  Figs.  5b  and  5c.  The  4  dB  pulse  width  is  44  nS  , 
and  the  total  pulse  compression  is  68  .  The  theoretical  compression  is 
3  x  23  or  69  .  The  sidelobes  are  not  as  small  as  theoretically  expected, 
nor  are  they  symmetric  about  the  central  peak.  This  is  due  to  the  fact  that 
(1)  the  variation  of  chirp  frequency  with  time  is  not  linear,  and  (2)  the 
amplitude  of  different  frequency  components  of  the  input  chirps  varies  by  as 
much  as  2  dB  (Fig.  5a). 

The  fact  that  the  output  pulse  width  is  not  limited  by  dispersion  is  a 
remarkable  result.  The  variation  in  group  delay  over  the  23  MHz  bandwidth 
is  .2  us  .  Thus,  one  might  expect  the  minimum  output  pulse  width  to  be  of 
this  order.  However,  it  is  shown  in  Appendix  A  that  dispersion  limits  Sezawa 
wave  convolvers  to  a  usable  time  bandwidth  product  of  roughly  100  .  In 


Appendix  A,  we  also  examine  the  effects  of  dispersion  on  correlation  and  show 
that  in  most  (but  not  all)  modes  of  operation,  dispersion  does  not  limit  the 
performance  of  the  storage  correlator. 

3.  Sezawa  Wave  Correlators 

A  series  of  storage  correlators  have  been  fabricated  and  used  in  storage 
applications;  however,  their  performance  is  inferior  to  the  best  of  our 
Rayleigh  wave  correlators.  A  number  of  these  problems  have  been  solved. 

The  first  problem  to  become  apparent  was  the  generation  of  bulk  waves  by 
our  IDTs.  This  is  more  of  a  problem  than  in  the  past  because  of  the  use  of 
4  finger  pair  (fp)  IDTs  rather  than  10  fp  IDTs.  Another  factor  is  that  the 
ZnO  films  are  thicker  (8  urn  for  Sezawa  wave  devices  and  1.6  urn  for 
Rayleigh  wave  devices),  and  so  bulk  waves  are  more  easily  excited  in  the 
ZnO/Si  structure. 

This  problem  has  been  solved  by  bonding  the  Si  substrate  to  a  brass 
carrier  box  with  indium  solder.  A  major  portion  of  the  bulk  wave  power  which 
reaches  the  back  surface  of  the  Si  substrate  is  transmitted  into  the  indium 
layer  and  brass  box  since  the  impedance  of  indium  and  silicon  are  about  the 
same  (ZSi  =  22,  Zjn  =  19)  ,  and  the  impedance  of  brass  is  slightly  higher 
(Zp  =  31)  .  Using  this  technique,  we  have  eliminated  the  problem  of  spurious 
bulk  waves  generated  by  the  interdigital  transdsucers. 

A  second  problem  is  associated  with  the  fact  that  our  top  plates  are  not 
rectangles  but  rhomboids  (Fig.  6a).  The  reason  for  using  slanted  ends  is  that 
most  of  the  spurious  surface  wave  power  generated  by  the  top  plate  writing  or 
reading  signal  misses  the  IDTs.  The  problem  with  using  this  design  for  Sezawa 
waves  is  that  the  SAW  velocity  is  significantly  different  in  the  metallized 
region  and  the  region  without  a  shorting  plane  at  the  surface.  If  one  is  only 


interested  in  delay  lines,  then  this  is  not  a  problem  if  both  ends  of  the  top 
plate  are  slanted  the  same  way  (Fig.  6a)  and  the  phase  distortion  at  one  end 
is  compensated  by  the  phase  distortion  at  the  other  end.  However,  if  the 
device  is  used  as  a  storage  correlator,  and  it  is  desired  to  store  charge  in 
diodes,  then  phase  distortion  of  more  than  45°  is  not  acceptable. 

For  Rayleigh  wave  devices,  this  has  not  been  a  problem  since  Av/v  is 
small  ( i . e . ,  the  SAW  velocity  in  the  shorted  region  is  about  the  same  as  the 
SAW  velocity  in  the  electrically  free  region).  For  a  25°  angle  on  the  top 
plate  and  1  mm  beamwidths,  the  phase  distortion  for  ZnO/Si  Rayleigh  wave 
devices  is  2°  .  For  our  ZnO/Si  Sezawa  wave  devices,  the  phase  distortion 
is  87°  !  Our  proposed  solution  is  shown  in  Fig.  6b.  This  eliminates  the 
phase  distortion  problems  and  also  should  decrease  our  electromagnetic 
feedthrough  problems  since  there  is  now  a  ground  plane  over  most  of  the 
device. 


C.  Monolithic  Storage  Correlator  Theory 

In  previous  progress  reports  and  papers,  we  have  included  a  great  deal  of 
experimental  data  on  the  dependence  of  correlator  output  on  various  parameters 
such  as  acoustic  voltage,  writing  signal  voltage,  etc.  We  have  developed  a 
comprehensive  theory  for  the  charging  and  reading  processes  in  the  monolithic 
storage  correlator.  A  complete  presentation  of  the  theory  is  given  in  a  paper 
accepted  for  publication  in  IEEE  Transactions  on  Sonics  and  Ultrasonics 
(Appendix  B).  Experimental  and  theoretical  results  are  given  for  the 
dependence  of  correlator  output  on  six  parameters:  acoustic  signal  duration 
and  amplitude,  write  signal  duration  and  amplitude,  and  read  signal  duration 
and  amplitude.  To  illustrate  this  theory,  we  present  two  of  these  results 


The  dependence  of  correlator  output  on  acoustic  signal  duration  is  shown 
in  Fig.  7.  Note  that  the  experimental  and  theoretical  results  are 
approximately  linear.  There  is  a  slight  upturn  in  both  the  experimental  and 
theoretical  results  when  the  acoustic  and  write  signals  have  approximately  the 
same  duration  (10  ms)  .  This  approximately  linear  characteristic  was 
surprising  when  it  was  first  experimentally  discovered  since  the  output  is 
very  nonlinear  when  both  the  write  and  acoustic  signals  are  varied  together. 
The  reason  for  the  linearity  is  clear  from  the  terms  in  the  theoretical 
expression.  When  the  write  and  acoustic  durations  are  varied  together,  the 
component  of  charge  with  no  spatial  variation  also  varies.  However,  if  the 
write  signal  duration  is  fixed  (as  in  Fig.  7),  then  the  component  of  charge 
with  no  spatial  variation  is  constant  (and  not  zero)  and  the  output  is 
proportional  to  the  component  of  charge  with  spacing  2tt/x  ,  i.e.,  to  the 
duration  of  the  acoustic  signal. 

The  dependence  of  correlator  output  on  write  signal  amplitude  is  shown  in 
Fig.  8  for  several  read  voltages.  The  output  saturates  for  write  voltages 
greater  than  1.5  V  (for  10  ms  charging  times). 

The  theory  has  been  used  in  the  design  of  improved  storage  correlators. 

It  has  also  been  used  to  predict  the  optimum  silicon  doping  density  to  use  for 
our  Sezawa  wave  devices  (Fig.  9).  The  optimum  doping  for  Sezawa  wave  devices 
(8  x  104cm-3)  is  twice  the  optimum  doping  for  Rayleigh  wave  devices 
(4  x  1014cm'3)  . 

D.  Comparison  of  Various  Systems  Based  on  Bandwidth  Considerations 

In  our  last  progress  report,  we  compared  the  theoretical  maximum 
bandwidth  obtainable  (when  the  optimum  tuning  circuit  is  used)  for  several 


different  systems  such  as  LiNb03  IDTs,  ZnO/Si  Sezawa  wave  IDTs,  and 
ZnO/Si  Rayleigh  wave  IDTs.  We  modelled  the  IDT  as  a  fixed  resistor  and 
capacitor  in  parallel  and  used  a  result  derived  by  Fano  to  compare  the 
systems. 

We  have  derived  a  theory  based  on  Fano's  work  to  predict  the  optimum 
bandwidth  for  an  IDT  when  the  IDT  is  represented  by  a  four  element  circuit 
model  instead  of  just  two  elements.  These  results  are  summarized  in  Appendix 
C.  The  differences  between  the  old  formulation  and  our  new  results  are  shown 
in  Fig.  10.  This  figure  is  a  plot  of  the  theoretical  maximum  bandwidth  as  a 
function  of  the  tuning  network  mismatch  loss.  As  can  be  seen,  the  inclusion 
of  a  frequency  dependent  radiation  resistance  and  reactance  dramatically 
limits  the  theoretically  maximum  bandwidth  (for  a  given  mismatch  loss). 

A  number  of  systems  of  interest  are  compared  in  Fig.  11.  LiNb03  IDTs 
are  still  the  best  system  to  use;  however,  it  can  be  seen  that  Sezawa  wave 
IDTs  are  almost  as  good.  A  number  of  additional  comparisons  are  made  in 
Appendix  C  such  as  the  dependence  of  bandwidth  on  the  number  of  finger  pairs, 
etc. 

E.  Schottky  Diode  Storage  Correlator 

In  the  past  six  months,  a  great  deal  of  progress  has  been  made  on  the 
Schottky  diode  storage  correlator.  We  have  demonstrated  such  a  correlator 
with  a  correlation  efficiency  of  -77.4  dBm  when  operating  in  the  acoustic  to 
plate  readout  mode.  More  importantly,  this  efficiency  was  achieved  using  a 
short  pulse  of  approximately  5  ns  duration  to  perform  the  storage 
operation.  Figure  12  shows  the  correlation  output  on  the  upper  trace  with  a 
one  volt  input  acoustic  signal,  a  one  volt  acoustic  readout  signal,  and 


63  dB  of  system  gain.  Note  the  second  correlation  output  resulting  from  the 
acoustic  readout  signal  reflecting  off  the  transducer  at  the  other  end  of  the 
device. 

The  lower  trace  shows  the  insertion  loss  of  the  system  for  a  one  volt 
amplitude  input  signal. 

Improvements  have  also  been  achieved  in  the  storage  times  of  these 

devices.  Whereas  our  earlier  devices  have  had  storage  times  of  the  order 

of  200  usee  ,  this  device  displayed  a  storage  time  close  to  1  msec  . 

Studies  we  have  carried  out  by  depositing  ZnO  on  large  area 
2 

(750,000  um  )  diodes  lead  us  to  believe  that  the  ZnO  sputtering  process 
degrades  the  storage  time  performance  of  the  diodes  by  a  factor  of  three  on 
the  average. 

The  largest  single  factor  responsible  for  the  increased  quality  of  these 
devices  is  the  improvement  of  the  PtSi  Schottky  barrier  technology.  We  have 
gained  access  to  an  extremely  clean  platinum  deposition  facility  owned  by 
Varian  Associates,  Inc.  This  is  capable  of  reducing,  quite  dramatically,  the 
level  of  trace  inpurities  at  the  Schottky  barriers,  which  would  tend  to 
increase  leakage  currents  and  lower  the  percentage  of  good  diodes  on  a 
particular  wafer.  We  have  indeed  observed  a  dramatic  improvement  in  the 
reliability  of  the  process  using  the  Varian  system.  The  next  step  to  be  taken 
in  improving  the  quality  of  these  devices  will  be  to  reduce  leakage  in  the 
diodes  due  to  the  high  electric  fields  at  the  edges  of  the  diodes,  the  so- 
called  "edge  effect"  leakage.  In  order  to  do  this,  we  propose  to  somewhat 
enlarge  the  diodes  so  that  it  will  be  possible  to  encircle  the  edge  of  each 
diode  with  a  p+  diffusion  "guard  ring,"  the  purpose  of  which  is  to  increase 
the  radius  of  curvature  at  the  edges  of  the  diode  and  hence  lower  the  electric 
fields  in  these  regions  and  hence  decrease  the  leakage  current  in  the  diodes. 


Storage  times  of  the  order  of  10  ms  are  expected  as  a  result. 

F.  Silicon  on  Lithium  Niobate 

We  have  also  studied  the  laser  annealing  of  polysilicon  deposited  on 
lithium  niobate  in  hopes  of  being  able  to  produce  semiconductor  quality  thin 
silicon  films  from  which  convolvers  and  correlators  can  be  made. 

Our  studies  have  been  carried  out  using  a  7  watt  argon  laser  with  a 
scan  system  capable  of  repeated  scans  of  a  40  pm  spot  over  a  large  area  of  a 
target. 

We  have  attempted  the  laser  annealing  on  lithium  niobate  samples  with 
barrier  layers  of  Si O2  ranging  in  thickness  from  .1  micron  to  1 
micron.  The  results  of  these  studies  indicate  that  it  is  indeed  possible  to 
anneal  the  small  grain  polysilicon  in  order  to  obtain  "long  grain"  polysilicon 
material  with  grain  dimensions  of  5  microns  wide  by  40  microns  long. 

These  grain  dimensions  would  prove  suitable  for  fabrication,  for  example,  of  a 
Rayleigh  wave  amplifier  as  a  monolithic  silicon  on  LiNbOj  structure. 

The  major  drawback  in  the  use  of  this  technique  at  the  present  time  is 
the  difficulty  in  preventing  the  generation  of  cracks  in  the  lithium 
niobate.  These  cracks  are  apparently  a  result  of  thermal  expansion  shock  in 
the  lithium  niobate  due  to  the  heat  absorbed  in  the  silicon  layer. 

The  situation  now  presented  to  us  is  one  of  finding  the  optimum 
conditions  of  laser  power,  Si02  layer  thickness,  and  silicon  deposition 
method  that  will  enable  us  to  anneal  the  silicon  without  damaging  the  lithium 
niobate. 


G.  FET  Storage  Correlator 

We  have  carried  out  theoretical  investigations  into  the  possibility  of 
making  a  device  similar  in  operation  to  the  storage  correlator  which,  like  the 
Schottky  diode  correlator,  will  have  the  advantage  of  fast  signal  read  in. 

The  basic  structure  of  this  device  is  shown  in  Fig.  13.  With  proper 
geometries  of  the  FET  array,  this  device  will  be  able  to  store  a  signal  in 
times  less  than  a  nanosecond,  which  is  quite  fast  enough  with  which  to  sample 
an  rf  cycle  that  is  8  ns  long. 

We  are  now  working  with  the  Integrated  Circuits  laboratory  to  design 


masks  and  a  processing  schedule  with  which  this  device  can  be  realized. 


Table  I 


Parameters  of  Sezawa  Wave  Convolvers 


ZnO  thickness 

7  pm 

SiO^  thickness 

2000  A 

Silicon  thickness 

350  pm 

IDT  and  top  plate  metal  thickness 

1000  A 

Beamwidth 

1  mm 

Length  of  interaction  region 

2.3  cm 

Silicon  orientation 

(100) 

IDT  finger  width  and  spacing 

8  pm 

Number  of  finger  pairs 

4 

Convolution  efficiency 

-50  dBm 

Insertion  loss 

31  dB 

Resonant  frequency 

155  MHz 

Bandwidth  (3  element  tuning  circuit) 

25  MHz 

Group  velocity  at  150  MHz 

4320  m/s 
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Figure  Captions 


1.  Dependence  of  insertion  loss  on  frequency  for  a  4  mm  long  Sezawa  wave 
delay  line.  The  5  element  IDT  tuning  circuit  is  also  shown. 

2.  Variation  of  convolution  efficiency  on  frequency  for  a  Sezawa  wave 
convolver.  The  3  element  IDT  tuning  circuit  is  also  shown. 

3.  (a)  5  bit  Barker  Code  correlation 
(b)  11  bit  Barker  Code  correlation 

4.  Circuit  diagram  for  chirp  compression  experiment. 

5.  Chirp  compression 

(a)  Input  chirp 

(b)  Output  pulse  (200  nS/di v) 

(c)  Output  pulse  (20  nS/di v) 

6.  (a)  Rayleigh  wave  top  plate  metallization 
(b)  Sezawa  wave  top  plate  metallization 

7.  Dependence  of  correlator  output  on  acoustic  signal  duration: 

_  Complete  theory;  _  Approximate  theory;  ....  Experimental 

8.  Dependence  of  correlator  output  on  write  signal  amplitude: 

_  Complete  theory;  _  Approximate  theory;  ....  Experimental 

9.  The  effect  of  doping  level  on  the  acoustic  wave  attenuation,  correlation 
efficiency  (without  attenuation)  and  total  correlation  efficiency 

10.  Comparison  of  bandwidth  predictions  using  the  full  transducer  model 

(  _  )  and  simple  R-C  model  (  _  )  for  a  Sezawa  wave  transducer  (fg  = 

165  MHz  ,  w  =  1  mm  ,  on  (100)  Si) 

11.  Comparison  of  UNbOj,  ZnO/Si  Sezawa  wave,  and  ZnO/Si  (first  Av/v 
peak  Rayleigh  wave  transducer  (fg  =  165  MHz  ,  w  *  1  mm) 
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12.  Schottky  diode  correlation  with  1  volt  input  signal  level 
Top  trace:  Correlation  output  200  mV/div 

Bottom  trace:  Insertion  loss  50  mV/div 

13.  A.  FET  correlator  structure  top  view 

B.  FET  correlator  structure  cross-sectional  view 
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APPENDIX  A 


A  Technique  for  Calculating  the  Bandwidth-Loss  Relationship  for  Acoustic  Transducers 


Ginzton  Laboratory 
Stanford  University 
Stanford,  CA  95305 


Abstract 

A  technique  is  developed  to  determine  the  maximum  bandwidth  ob¬ 
tainable  from  a  transducer  matched  through  the  optimum  matching  network 
to  an  ideal  source.  The  transducer  is  modeled  by  four  elements.  This  tech¬ 
nique  can  be  applied  to  surface  acoustic  wave  (SAW)  transducers  or  bulk 
acoustic  wave  transducers.  The  technique  was  used  1)  to  compare  SAW  IDTs 
on  various  materials  such  as  LiNb03  and  ZnO/Si,  2)  to  compare  different 
modes  such  as  the  Rayleigh  and  Sezawa  modes  in  ZnO/Si  structures,  and 
3)  to  compare  different  substrate  orientations.  It  is  found  that  Sezawa  wave 
IDTs  in  ZnO/(100)Si  offers  the  best  performance  of  the  monolithic  structures 
considered,  and  that  these  transducers  are  as  good  as  LiNb03  IDTs  for  ob¬ 
taining  broad  bandwidth.  The  design  of  transducers  for  maximum  bandwidth 
is  also  considered  using  this  technique. 


L  Introduction 


A  major  problem  in  the  design  of  SAW  devices  is  obtaining  the  broad¬ 
band  operation  necessary  for  many  applications,  such  as  utilization  in  spread 
spectrum  systems.  There  has  been  a  major  effort  in  developing  monolithic 
versions  of  many  SAW  devices  such  as  delay  lines,  convolvers,  and  storage 
correlators.  These  devices  have  been  narrowband,  and  consequently,  there 
was  a  search  for  other  monolithic  systems  which  had  broader  bandwidth 
such  as  Sezawa  wave  IDTs  and  edge  bonded  transducers.  These  systems  have 
been  compared  by  considering  the  SAW  coupling  coefficient  or  transducer  Q; 
however,  the  transducer  bandwidth  (for  a  given  loss)  is  a  function  of  several 
of  these  parameters,  and  it  is  not  always  evident  which  system  is  best.  In  this 
paper,  a  technique  is  developed  to  theoretically  compare  such  alterntives  on 
the  basis  of  bandwidth  considerations.  After  a  choice  of  transducer  is  made, 
the  technique  described  here  can  be  used  to  determine  the  optimum  design 
criteria.  Smith  et.  al.  [1]  and  Khuri-Yakub[2]  have  addressed  the  problem 
of  transducer  design  and  hcuristically  argued  that  the  acoustic  Q  should  be 
equal  to  the  electronic  Q  for  optimum  bandwidth.  This  is  an  approximate 
result  which  is  only  valid  for  small  mismatch  losses  and  series  inductor  tuning. 
This  rule  is  not  applicable  to  more  sophisticated  tuning  circuits.  The  tech¬ 
nique  described  in  this  section  determines  the  optimum  transducer  parameters 
assuming  the  best  matching  network  for  these  parameters  is  used.  In  par¬ 
ticular,  balanced  and  unbalanced  operation  can  be  compared,  and  the  op¬ 
timum  width  and  number  of  finger  pairs  can  be  determined. 


EL  Theory 


1.  History 


Bode  [3]  recognized  that  a  source  can  be  matched  to  a  load  consisting  of 
a  capacitor  and  a  resistor  in  parallel  (Fig.  1)  only  over  a  limited  frequency 
range.  He  obtained  the  limitation  on  the  bandwidth 


|p(w)irfW  ^  RC 


where  p  is  the  input  reliection  coefficient. 


p(oj)  = 


Z(u)  -  1 
Z(u)  +  1 


(1) 

(2) 


where  Z(tu)  is  the  impedance  of  the  matching  network  and  load.  If  p  is 
constant  (pm)  over  the  desired  bandwidth^#)  and  equal  to  1  over  the  rest  of 
the  frequency  spectrum,  then  Eq.  (1)  becomes 
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As  expected,  the  price  of  broader  bandwidth  is  higher  loss. 

As  a  first  approximation  the  circuit  model  of  an  IDT  can  be  represented 
as  a  fixed  value  resistor  and  capacitor  in  parallel,  and  Eq.  (1)  can  then  be 
used  to  determine  the  relationship  between  bandwidth  and  loss.  This  has 
been  done  by  Reeder[4]  and  Schellenberg[5j. 

The  problems  with  this  approximation  (see  Fig.  2a)  are  1)  the  acoustic 
conductance  is  not  at  all  constant,  and  2)  the  acoustic  susceptance  is  nonzero: 


CM  = 


X 


B(u) 


f,  =  Gsjn(2X)-2X 


2X2 


(4) 

(5) 

1 


where 


(6) 


x  _  N ?r(w  —  wQ) 

Wo 

and  Wo  is  the  resonant  frequency  and  N  is  the  number  of  finger  pairs. 

Fano[6j  extended  the  work  of  Bode  by  considering  the  matching  of  a 
source  to  an  arbitary  passive  load.  He  then  applied  his  results  to  the  case  of 
a  resistor  and  inductor  in  parallel  (the  dual  network  of  Bode’s  network)  and 
rederived  Eq.  (1). 

In  the  next  section,  we  shall  demonstrate  that  a  passive  circuit  composed 
of  four  circuit  elements  with  frequency  independent  values  can  represent  the 
IDT  admittance  represented  by  Fig.  2a  and  then  use  this  model  with  Fano’s 
theory  to  obtain  the  bandwidth-loss  relationship  in  the  general  case  for  an 
acoustic  transducer. 

2.  Circuit  Model 

The  circuit  model  shown  in  Fig.  2b  has  been  used  to  model  bulk  wave 
transducers  [7], -however,  this  is  the  first  instance  the  author  is  aware  of  that 
the  model  has  been  applied  to  SAW  transducers. 

The  relation  between  the  parallel  circuit  for  the  IDT  (Fig.  2a)  and  the 
model  (Fig.  2b)  can  be  found  by  expanding  the  expressions  for  the  transducer 
admittance  Eq.  (4)  and  (5)  about  the  transducer  center  frequency,  and  ex¬ 
panding  the  model  admittance  about  the  same  frequency.  If  the  zero  and  first 
order  terms  in  (w  — wo)  for  the  real  and  the  imaginary  parts  of  the  transducer 
admittance  are  equated,  then  there  are  four  equations  to  determine  the  four 
model  parameters:  R,Cq,L  and  Cx 

S  =  on 

t Jao 

Cq  ~  Cj  -f-  Cp  (8) 
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(9) 


Ni r  _  Cp 
SuJoGao  2G^o 

Cl  =  LuJ 

where  CT  is  the  transducer  finger  capacity  and  Cp  is  the  pad  and  stray 
capacity.  For  most  SAW  transducers,  Av/v  is  small  enough  (Av/v  <  .03) 
that  the  imaginary  part  of  the  admittance  is  predominantly  due  to  the  parallel 
capacitor,  and  a  better  match  between  the  transducer  and  model  admittances 
is  obtained  if  the  first  three  terms  in  the  expansion  of  the  real  part  of  the 
admittance  are  equated,  and  only  the  first  term  of  the  imaginary  parts  is 
equated.  The  model  relations  in  this  case  are  the  same  as  Eqs.  (7)  to  (10) 
except  that  the  expression  for  L  becomes 

T  _  nN 

L  —  - 

2V3ci?o  Gao 

For  the  case  of  a  ZnO/Si  Sezawa  ware  transducer  the  real  and  imaginary  parts 
of  the  transducer  admittance  are  shown  in  Fig.  3.  Also  shown  in  these  figures 
are  the  model  predictions.  Note  that  the  agreement  is  good  near  the  resonant 
frequency,  but  for  larger  bandwidths,  the  real  part  of  the  model  admittance 
does  not  drop  as  rapidly  as  the  transducer  admittance.  Hence,  the  predictions 
made  in  the  next  section  will  be  optimistic  for  large  bandwidths.  This  was  the 
case  for  all  of  the  types  of  transducers  studied.  For  the  sake  of  comparison, 
the  transducer  and  model  parameters  were  determined  for  1)  Sezawa  wave 
IDTs:  ZnO/  (100)Si  and  (111)  Si,  2)  First  order  Rayleigh  wave  IDTs:  ZnO/ 
(100)Si  and  (lll)Si  for  the  first  and  second  peaks  in  Av/v  and  3)  Rayleigh 
wave  IDTs  on  YZ  LiNb03.  These  parameters  are  tabulated  in  Table  1  for 
use  in  a  later  section  of  this  paper.  The  capacity  of  the  bonding  pads  Cp 
is  neglible  for  transducers  on  a  semi-infinite  LiNb03  substrate.  However,  all 
of  the  ZnO/Si  monolithic  transducers  considered  here  have  either  IDTs  at 
the  top  surface  and  a  shorting  plane  at  the  interface  (Sezawa  wave,  Rayleigh 
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wave  (first  Av/v  peak))  or  IDTs  at  the  interface  and  a  shorting  plane  at  the 
top  surface  (Rayleigh  wave  (second  Av/v  peak)).  Consequently,  the  bonding 
pad  capacity  cannot  be  neglected  for  monoiithic  devices. 

3.  Integral  Equations 

Consider  the  system  in  Fig.  4a.  Any  physically  realizable  load  can 
be  represented  by  a  reactive  two  terminal  network  terminated  by  a  resistance. 
This  resistance  can  be  made  equal  to  the  source  resistance  by  incorporating  a 
transformer  in  the  reactive  network  (Fig.  4b).  It  is  obvious  from  Fig.  4b  that 
the  transducer  network  will  limit  the  power  transmission  from  the  source  to 
tne  load  resistor;  independent  of  what  the  matching  network  (Nm)  is.  Fano’s 
analysis  is  more  tractable,  if  the  two  networks  are  flipped  (Fig.  4c)  so  the 
source  is  driving  the  transducer  network  Ny.  Fano  considers  the  function 
lnl /p  where  p  is  the  reflection  coeffient  of  the  matching  network  and  load, 
and  shows  that  for  each  zero  in  the  right  half  plane  of  the  load  transmission 
coefficient  ti,  an  integral  equation  i3  obtained  relating 

•  p  (reflection  coefficient  of  the  load  plus  matching  network) 

•  u  (angular  frequency) 

•  X,  (zeros  of  the  matching  network  that  lie  in  the  right  half  plane.) 

In  particular,  he  obtained  the  relations 

.r^2“+!:tasd“=(-i)i^+'  <i2) 

for  a  zero  of  order  k  at  the  origin,  and 

=  (13) 

for  a  zero  of  order  k  at  infinity.  A*  is  the  ith  coefficient  of  the  Taylor  series 
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expansion  of  lnl/|pi|  about  the  point  j  and 


^i+1  =  ^+,  -  Jj^TJ  E  Xr(”+,)  (14) 

F&+ 1  =  ^£+1  -  3+7  E  >f+I  (15) 

'  t 

The  preformance  of  the  system  is  given  by  the  left  hand  side  of  the 
integral  equations  and  the  limitations  on  the  performance  is  given  by  the 
right  hand  side.  The  transducer  characteristics  enter  through  Aj  and  the 
matching  network  characteristics  enter  through  the  positions  of  the  zeros  of 
the  matching  network  X,-. 


4.  Application  to  Acoustic  Transducers 


The  impedance  of  the  model  shown  in  Fig.  2b  is 

7  _  11  +  s(Ci  +-  Co)  +-  s>‘LCx  +•  $*LCxCo 
T  ~  R  sCi  +-  s2C\Co 


(16) 


where  s  =  ju.  All  impedances  have  been  normalized  to  Rs.  The  reflection 
coefficient  is 


Pi  = 


Zt  —  1 
Zt  +■  1 


(17) 


The  transmission  coefficient  is 


= _ 2(sRC1-\-s2RC,Cq) _ 

1  1  +  s(Cx +- C0  +  flC!)  +  s2(L(7i  +  RCxCo)  +  s*LCiC0  {  ’ 

There  is  one  zero  of  t\  at  the  origin  and  two  zeros  at  infinity.  The  function 
of  interest  is 

_L  =  I  1  +  «(Ci  +-  Co  +  RCi)  +-  s2(LCx  +  RCiCq  +■  s^LCjCo 
l/jil  ni  +  s(C1  +  Co-HC1)  +  s2(LC1-i?C16,o)  +  s3LC1C0 


Expanding  F  about  the  origin  yields 


F  —  Ajt*,'  -|- 


(20) 


Thus, 


A?  =  ~|w==0  =  2  RCi 
dw 


Expanding  F  about  infinity  is  equivalent  to  expanding  F  about  q 

q-0. 

F=Af9+Afy +  ... 


(21) 
1  /u  as 

(22) 


Thus, 


2  R 
L 


where 


1  dzF . 

6d^l?=0  “ 


2 R  .  _  R2C3 
L2C3{  3 L  } 


1 

C3 


Note  that  C3**Ci  and  from  Eqs.  (7)  to  (10), 


£i  1 

L  ™{NR)2 


(23) 

(24) 


(25) 


(26) 


so 

A~—  ‘27> 

For  an  IDT,  the  number  of  finger  pairs  N  is  at  least  one;  hence,  A£°  is  negative. 
This  is  in  general  true,  and  as  we  shall  see,  this  is  important  for  a  solution 
to  these  equations  to  exist. 

Using  Eqs.  (12), (13), (21), (23)  and  (27)  we  obtain  the  results 


;?lnj^  =  f(2RC‘-2?*rI) 


(28) 


8 


(29) 


roo 

Jo 


.  1  ,  “K .  2  R  v—v  . 

In  — du>  =  — ( - 2  >  S/) 

|p|  2  KL  *} 


f 


2,  1  ,  7T .  3  R  .  1  .  2  . 

“  ,nR‘,“=2lI5^(1_^)+^i‘1 


(SO) 


Examine  first  Eq.  (29).  This  equation  is  the  dual  of  the  equation  obtained 
by  Bode  when  $,•  =  0.  Good  matching  occurs  when  p**0  and  lnl/|p|  is 
large.  Thus,  the  left  hand  side  of  Eq.  (29)  is  the  area  under  a  matching 
versus  frequency  plot.  The  value  for  the  area  depends  on  the  transducer 
parameters  and  will  be  largest  when  s*  =  0.  However,  since  Eq.  (28)  and 
Eq.  (30)  must  also  be  satisfied,  there  will  have  to  be  several  nonzero  roots  st-. 
For  the  applications  of  interest  in  this  paper,  a  fiat  frequency  response  over 
a  given  bandwidth  is  desired.  According  to  Eq.  (29),  the  best  transmission 
over  this  bandwidth  will  occur  when  the  transmission  outside  the  bandwidth 
is  zero  and  no  "area”  is  wasted  outside  the  desired  bandwidth.  Thus  the  ideal 
frequency  dependence  of  the  reflection  coefficient  is 


%K  whenwi  <  w  <  w2 
1  otherwise 


Eqs.  (28)  to  (30)  become 


(31) 


f(—  )  =  RCi--~- 

2  Wi  W2  52 


(32) 


K.  R 

y  (w2  —  Wl)  =  —  —  Si  —  S2  —  .  .  . 


K .  3  3.  3R 

T«— 


£  +  4+4 


(33) 

(34) 


For  the  matching  network  to  be  realizable,  the  zeros  in  the  right  half  plane 
will  be  symmetric  about  the  real  axis.  Thus  the  st-  are  real  or  occur  in 
complex  conjugate  pairs.  We  shall  assume  that  the  transducer  parameters 
(/?,  L,  Ci,  Co)  and  the  desired  passband  frequencies  (o/i  to  a>2)  are  known 
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and  that  it  is  desired  to  find  the  matching  parameter  K.  At  least  three 
variables  if,si,s2  are  needed  to  solve  the  three  equations.  Are  more  zeros 
needed  for  the  optimum  match?  If  N  additional  zeros  are  included  and  N 
additional  minimizing  equations  are  used;  then  it  is  easily  shown  that  K  is 
a  monatonically  decreasing  function  of  N.  In  other  words,  the  best  match 
occurs  when  only  enough  zeros  are  included  to  simultaneously  solve  the  three 
equations  (34)  to  (36).  In  our  case  two  zeros  are  needed. 

If  we  recall  that  the  real  parts  of  s*  are  positive  and  that  these  zeros  are 
either  both  real  or  complex  conjugate  pairs,  then  it  can  be  seen  that  Eq.  (32) 
and  Eq.  (33)  set  maximum  values  for  K 


2  R 


K 


i  _ 


L(w  2  —  <*>l) 
2jRCiWlW2 


U>2  —  k'l 

Similarly,  Eq.  (34)  provides  the  minimum  for  K 


when  Si  =  s2  =  0 


when  si  =  s2  =  00 


(35) 

(36) 


K-min  — 


6  R  _  2R 
LaCi  17 

«i-«5 


(37) 


The  value  of  K  which  satisfies  all  three  equations  is  between  these  limits. 
Using  Eq.  (7)  to  Eq.  (10),  it  can  be  seen  that  K*max  is  slightly  smaller  than 


*maz 


TC1  z? 

**■  m/rv  _  _  /  D  v  O  -  . 

=  LCiWiW2^1  —  (— )  —  1 


v  ,  (38) 

where  B  is  the  fractional  bandwidth.  A  solution  for  K  can  only  exist  if 
Krnin  ^ 

£sl1  —  1  d2  r>.  cn3 

i 

Thus, 

(40) 


,  ^  IC'na,  _  g(l  -  £)  1  B*w  <7,  6B* 

^K—  -  1  -  zn  ~(1  +  3^  ~  T))fl  ~  CbM  ~  W 


(39) 


B  <  — 
N 


1 
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This  is  equivalent  to  the  intuitive  argument  that 


B  < 


_1_ 

N 


(41) 


i.e.  the  bandwidth  is  approximately  the  frequency  range  between  the  3dB 
points  of  the  conductance  of  the  transducer.  As  might  be  expected,  the 
maximum  bandwidth  is  larger  than  this  value  by  \/6  times  this  value.  A 
simultaneous  solution  to  the  three  Eqs.  (32)  to  (34)  was  obtained  using  a 
computer.  The  results  are  given  in  the  next  section  for  a  number  of  different 
situations. 

Three  comments  regarding  bulk  acoustic  wave  transducers  are  in  order. 
First,  the  expressions  for  the  model  parameters  are  different.  If  Rb  + 
1/joJoCs  is  the  transducer  impedance  at  center  frequency,  then  the  model 
parameters  are  [8,9] 


R  =  1 

Rb(uoCb)2 

(42) 

Cq  =  Cb 

(43) 

k2ff 

<?i  =  ~ — 

1  -  k*ss 

(44) 

L  ~  Wc<;i 

(45) 

?*- 
«  tO 

s 

II 

nJS 

(46) 

where 


Secondly,  unlike  the  result  with  SAW  tranducers,  it  is  found  that  for  bulk 
transducers  i K'max  >  Kma z.  Finally,  in  some  cases  a  flat  frequency  response 
is  not  important;  rather,  a  good  pulse  response  is  desired.  Thus,  a  gaussian 
frequency  response  is  a  better  choice,  and  a  relation  between  the  minimum 
loss  and  the  width  of  the  gaussian  is  desired.  An  interesting  result  is  obtained 
in  this  case,  namely  that  the  left  hand  side  of  two  of  the  integrals  diverge.  In 


other  words,  it  is  not  possible  to  match  a  transducer  for  a  gaussian  bandpass 
over  either  of  the  ranges  Wi  to  oo  or  u>2  to  0  (dc)  where  Wi  <  ujq  <  w2-  If  a 
finite  cutoff  is  assumed  on  both  ends  of  the  response,  then  two  of  the  integrals 
can  be  carried  out  exactly,  and  a  good  approximate  closed  form,  expression 
can  be  obtained  for  the  third  integral.  The  analysis  for  bulk  acoustic  wave 
transducers  can  then  be  carried  out  just  as  for  SAW  IDTs. 


CL  Results  and  Discussion 

Has  all  of  this  formalism  done  any  good,  or  should  we  just  have  used 
Bode’s  result?  The  difference  in  Bode’s  result  and  our  result  is  shown  in  Fig. 
5  for  Rayleigh  and  Sezawa  wave  transducers.  Obviously,  Bode’s  result  is 
not  very  applicable  to  SAW  transducers.  A  number  of  different  systems  and 
transducer  designs  will  now  be  compared. 

Interdigital  transducers  using  1)  1st  Rayleigh  mode  (hereafter  referred 
to  as  Rayleigh  wave)  on  YZ  LiNb03  ,  2)  Rayleigh  waves  on  ZnO/(100)Si 
and  (lll)Si,  and  3)  Sezawa  waves  on  ZnO/(100)Si  are  compared  in  Fig. 
6.  ZnO/Si  Sezawa  wave  IDTs  are  theoretically  slightly  superior  to  LiNb03 
IDTs  due  to  the  higher  value  for  Av/t’  (.023  versus  .024).  In  practice,  the 
full  bandwidth  tor  Sezawa  wave  IDTs  may  be  difficult  to  achieve  due  to 
the  significant  problems  with  stray  capacity.  It  is  important  for  broadband 
matching  to  any  transducer  to  keep  the  stray  capacity  much  smaller  than 
the  IDT  series  capacity.  The  series  capacity  of  ZnO/Si  Sezawa  wave  IDTs 
is  much  smaller  than  LiNb03  IDTs  (.45  pF  versus  1.3  pF  for  the  geometry 
considered  here)  since  the  dielectric  constant  of  ZnO  is  five  times  smaller  than 
LiNb03.  Consequenty,  .5  pF  stray  capacity  in  the  connecting  leads  is  much 
more  detrimental  for  ZnO/Si  IDTs.  This  result  is  shown  in  Fig.  6,  where  it 
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can  be  seen  that  LiNb03  IDTs  are  superior  for  mismatch  losses  greater  than 
ldB. 


Another  important  prediction  is  that  the  use  of  (100)Si  instead  of  (lll)Si 
for  Rayleigh  -wave  monolithic  devices  should  produce  broader  bandwidth 
devices.  Traditionally,  (lll)Si  has  been  used  for  monolithic  SAW  devices 
at  Stanford;  however,  as  these  results  show,  broader  bandwidth  devices  can 
be  obtained  on  (100)Si.  The  physical  reasons  for  the  superiority  of  (JOO)Si 
are  1)  Av/v  is  larger  (.0046  versus  .0039  )  and  2)  the  peak  in  Av/v  occurs  at 
a  smaller  value  of  hk  (.31  versus  .37)  and  consequently,  the  series  capacity  is 
larger,  and  stray  capacities  are  less  significant. 

Six  kinds  of  monolithic  transducers  are  compared  in  Fig.  7.  Sezawa  wave 
transducers  are  slightly  better  than  second  peak  Rayleigh  wave  transducers. 

In  practice,  the  performance  of  the  second  Av(v  peak  Rayleigh  wave  transducers 
will  be  significantly  degraded  since  the  transducer  is  at  the  ZnO/Si02  in¬ 
terface,  and  the  calculation  for  the  series  and  pad  capacities  assumed  an 
infinately  thick  Si02  layer.  Actually,  the  layer  would  only  be  several  microns 
thick,  and  a  larger  parallel  capacity  would  be  present  due  to  conduction  in  the 
silicon.  Note  that  this  is  parallel  capacity  (shorting  the  radiation  resis¬ 
tance),  not  series  capacity  since  thi3  capacity  is  due  to  fields  between  fingers 
which  dees  not  generate  acoustic  stress  since  Si02  and  Si  are  not  piezoelectric. 

Balanced  and  unbalanced  operation  is  compared  in  Fig.  8.  As  expected 
from  calculations  of  electronic  Q,  balanced  operation  is  better.  However,  it  is 
interesting  to  note  that  unbalanced  operation  is  better  for  large  bandwidths 
(and  high  loss).  Also,  the  increase  in  bandwidth  when  switching  from  un¬ 
balanced  to  balanced  operation  i3  not  2:1  as  expected  from  the  increase  of  a 
factor  of  2  in  electronic  Q.  This  is  well  verified  experimentally. 

The  physical  reason  for  the  results  listed  in  the  above  paragraph  is  the 
tradeoff  between  1)  electronic  limitations  to  the  bandwidth  due  to  the  size  of 


the  finger  capacity,  and  2)  acoustic  limitations.  As  a  rough  approximation^], 
we  can  equate  the  acoustic  and  electronic  Q’s  to  determine  the  optimum 
number  of  finger  pairs: 

Qa  =  N  (47) 


Qe  =  — L_  =  — - — L 

wCtRo.0  2.13^2-  N 


(48) 


where  v  =  1  for  balanced  operation,  and  v  =  2  for  unbalanced  operation. 
Equating  Q’s 


N  =  (- 


v 


(49) 


2.13^ 

For  ZnO/Si  Sezawa  ware  balanced  transducers,  the  optimum  N  according 
to  this  argument  is  4  while  for  LiNbC>3  IDTs  the  optimum  N  is  4.5.  The 
prediction  of  the  theory  based  on  Fano’s  work  is  shown  in  Fig.  9.  It  can  be 
seen  that  the  optimum  N  is  3.5  for  both  ZnO/Si  Sezawa  wave  and  LiNbOj 
IDTs. 


IV.  Conclusions 

A  powerful  and  versatile  technique  has  been  developed  in  this  paper  for 
comparing  different  transducer  systems.  Sixty  per  cent  bandwidth  with  low 
mismatch  loss  (<3dB)  is  predicted  for  Sezawa  wave  IDTs.  This  technique 
can  be  used  for  transducer  design,  and  the  optimum  number  of  finger  pairs, 
optimum  width,  and  optimum  finger  width  to  spacing  can  be  determined. 
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Table  1.  Parameters  of  Systems  Considered  in  this  Paper 


System 

Mode 

Gao 

Ct 

cP 

R 

Co 

Ci 

(mmho) 

(PF) 

(PF) 

(n) 

(pF) 

(pF) 

ZnO/(100)Si 

Sezawa 

.145 

.401 

.1 

6890. 

.501 

.0386 

ZnO/(lll)Si 

Sezawa 

.114 

.388 

.08 

8770. 

.468 

.0303 

ZnO/(100)Si 

Rayleigh  * 

.0834 

1.40 

.42 

12000. 

1.82 

.0222 

ZnO/(lll)Si 

Rayleigh  * 

.0591 

1.16 

.35 

16900. 

1.51 

.0157 

ZnO/(100)Si 

Rayleigh  -f~ 

.0758 

.485 

.063 

13200. 

.548 

.0201 

ZnO/(lll)Si 

Rayleigh  + 

.0954 

.478 

.065 

10500. 

.543 

.0253 

YZ  LiNbOa 

Rayleigh 

.222 

1.29 

0. 

4500. 

1.29 

.0787 

*  First  A*  peak 
-f  Second  ^  peak 
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THEORETICAL  AND  EXPERIMENTAL  RESULTS  FOR  MONOLITHIC  SAW  MEMORY  CORRELATORS 
H.  C.  Tuan,  J.  E.  Bowers,  and  G.  S.  Kino 


ABSTRACT 

A  detailed  theory  is  presented  for  the  writing  and  reading  processes  in 
the  monolithic  surface  acoustic  wave  (SAW)  storage  correlator.  Theoretical 
and  experimental  results  are  presented  on  the  dependence  of  the  correlation 
output  on  the  duration  of  the  acoustic,  write,  and  read  signals.  The 
dependence  of  the  output  on  the  amplitudes  of  the  acoustic,  write,  and  read 
signals  is  also  given.  The  theoretical  predictions  are  in  good  agreement  with 
the  experimental  results  in  all  of  these  cases.  Suggestions  are  provided  for 
the  design  of  improved  monolithic  storage  correlators. 


I.  Introduction 

The  diode  storage  correlator  was  first  implemented  using  an  air-gap 
structure  where  the  surfaces  of  a  LiNb03  and  a  Si  substrate  were  placed  a  few 
thousand  Angstroms  apart. A  comprehensive  theoretical  and  experimental 
study  of  this  device  was  undertaken  by  Borden  and  Kino®*^  and  Loh  et  al.5 
El -Nokali  and  Adler®  extended  the  portion  of  the  theory  dealing  with  the 
charging  process  to  apply  to  diodes  with  minority  carrier  lifetimes  longer 
than  the  rf  period. 

The  air-gap  correlator  has  a  number  of  problems  such  as  spurious  bulk 
wave  generation,  air-gap  uniformity,  and  the  necessity  for  large  voltages  to 
store  and  read  out  signals.  Consequently,  Tuan,  Khuri  -Yakub,  and  Kino?*® 
developed  the  monolithic  storage  correlator  which  used  the  first  order 
Rayleigh  mode  propagating  in  a  ZnO-on-Si  structure.  This  reduced  or 
eliminated  the  problems  described  above.  However,  a  limitation  of  this 
approach  is  that  narrower  bandwidths  are  obtained  than  with  LiNbOj 
devices.  Monolithic  devices  with  broader  bandwidths  may  be  obtained  through 
the  use  of  the  second  Rayleigh  mode.^-H  The  results  presented  in  this  paper 
are  from  a  device  which  used  the  first  Rayleigh  mode. 

A  complete  theory  of  the  monolithic  storage  correlator  is  presented 
here.  The  results  of  this  theory  are  compared  to  a  comprehensive  set  of  data 
on  the  dependence  of  the  correlator  output  on  six  parameters:  acoustic  signal 
duration  and  amplitude,  write  signal  duration  and  amplitude,  and  read  signal 
duration  and  amplitude.  Results  are  presented  for  both  the  normal  correlation 
mode,  when  the  acoustic  and  plate  signal  durations  are  of  the  order  of  the  top 
plate  duration  (typically  several  microseconds),  and  for  the  input  correlation 
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mode  where  the  acoustic  and  plate  signal  durations  are  of  the  order  of  the 
storage  time  of  the  diode  (typically  hundreds  of  milliseconds).  Only  the 
correlation  of  two  square  pulses  is  considered  here.  Other  applications  such 
as  the  use  of  chirp  or  phase  shift  keyed^^  signals  involve  simple  calculations 
of  compression  gains. 

II.  Description  of  SAW  Memory  Correlator 

A  schematic  drawing  of  the  monolithic  storage  correlator  is  shown  in 
Fig.  1.  The  operation  of  the  correlator  is  a  two-step  operation  where  the 
correlation  Q(z)  of  the  acoustic  signal  VA(t)  from  the  left  transducer  and 
the  writing  signal  Vy(t)  are  first  stored'  in  the  diode  array  in  the  form  of 
a  spatially-varying  charge  pattern.  Then  a  second  reading  signal  VR(t )  is 
applied  to  the  center  (or  plate)  port.  As  a  result  of  the  interaction  of  the 
stored  charge  pattern  and  the  plate  signal,  the  correlation  and  convolution  of 
the  stored  signal  Q(z)  and  the  read  signal  V^(t)  are  obtained  at  the  left 
and  right  interdigital  transducers  (IDT's),  respectively.  A  second  method  of 
reading,  called  acoustic -to -pi ate  reading,  is  a  simple  extension  of  the  theory 
presented  below. 

The  parameters  of  the  device  used  in  the  experiments  described  here  are 
given  in  Table  I.  The  parameters  in  the  first  part  of  the  table  are  measured 
parameters,  and  the  second  part  of  the  table  contains  calculated  parameters 
sych  as  the  equivalent  acoustic  capacitance.^  The  definitions  of  the 
variables  used  in  this  paper  are  also  given  in  Table  I.  In  particular,  the 
capacity  of  the  diode  at  zero  bias  is  C<j  ,  the  capacity  of  the  region  in 
between  the  diodes  is  ,  the  total  silicon  capacity  is  CQ  »  +  C<j'  , 
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the  capacity  of  the  ZnO  and  SiC^  layers  is  Cox  ,  and  the  total  capacity  of 
the  plate  is  Cp  -  C0xC0/(C0x  +  C0)  •  A11  of  the  capacitances  and  charges  are 
per  unit  area.  Capital  subscripts  are  used  to  designate  external  voltages 
such  as  ,  and  the  internal  voltage  is  given  by  lower  case  subscripts 

(Va)  • 

The  equivalent  circuit  model  used  in  the  writing  process  is  given  in 
Figs.  2a  and  2b.  This  model  has  been  changed  from  that  used  earlier  for  the 
air-gap  correlator^ 6  for  two  reasons.  The  substrate  resistance  cannot  be 
neglected  for  devices  which  use  silicon  with  resistivities  of  the  order  of 
10  Ocm  ,  unless  epitaxial  layers  on  highly  conducting  substrates  are  used. 
Secondly,  it  is  important  to  realize  that  the  acoustic  wave  does  not  interact 
directly  with  the  plate  circuit  so  that  neither  part  of  the  circuit  loads  the 
other.  It  is,  therefore,  more  correct  to  use  two  separate  circuits  for  the 
two  branches,  as  shown  in  Fig.  2a  and  2b. 

Because  of  the  high  density  of  diodes,  the  resistance  of  the  plate 
circuit  can  be  calculated  without  including  the  effects  of  spreading 
resistance. 

When  the  dc  bias  on  the  plate  is  varied,  the  MOS  capacity  of  the  region 
in  between  the  diodes  C^'  also  varies,  and  the  rf  coupling  of  the  plate 
and  acoustic  ports  to  the  diode  changes.  Due  to  charge  injection  in  the  ZnO, 
the  dependence  of  the  MOS  capacity  on  the  dc  bias  is  complicated.  This 
effect  has  been  studied  by  Pierret  et  al.1A 

The  model  shown  in  Fig.  2a  is  for  a  single  diode,  and  the  source  voltage 

vp  is  the  voltage  at  the  top  plate.  The  external  driving  voltage  Vw  is 
connected  to  the  top  plate  through  a  short  load.  The  rf  impedance  of  a 
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short  lead  (~  1  cm)  can  easily  be  of  the  same  order  of  magnitude  as  the 

impedance  of  the  plate  capacitor,  and  thus  this  inductance  L  cannot  be 
neglected. 

The  values  of  the  electrical  half  of  the  circuit  (R,  L,  and  Cp)  can  be 
determined  by  measuring  the  plate  impedance  as  a  function  of  frequency 
(Fig.  3).  The  impedance  of  the  circuit  in  Fig.  2a  matches  the  measured 
impedance  quite  well.  The  values  used  in  the  fit  (see  Fig.  3)  agree  with  the 
component  values  expected  from  calculations  of  R,  L,  and  Cp  (see  Table  I). 

The  circuit  models  used  for  the  diode  (Fig.  2b)  and  the  read-out  process 
(Fig.  2c)  are  the  same  as  the  models  used  for  the  air-gap  correlator. ^  in  the 
read-out  process,  the  capacitance  of  an  individual  diode  is  regarded  as  being 
changed  by  the  writing  process.  When  an  rf  signal  is  applied  to  the  platq, 
it  develops  a  potential  across  this  reverse  biased  diode  capacity,  which 
depends  on  its  capacity  and  on  the  capacity  in  series  with  the  diode. 

III.  Theory 

We  shall  initially  consider  the  writing  process  and  calculate  the 
component  of  stored  charge  with  wave  number  k  .  Then,  the  read-out  process 
is  considered,  and  the  output  voltage  is  calculated. 

Writing 

The  model  shown  in  Fig.  2a  is  essentially  just  a  sample  and  hold  circuit, 
and  the  effect  of  the  resistances,  inductance,  and  parallel  capacitance  C^j 
is  to  reduce  the  voltage  at  the  diode  V<j  . 
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In  the  first  few  cycles  of  operation,  when-  the  rf  potential  becomes 
positive,  the  diode  potential  rises  to  approximately  0.5V  .  During  the 
negative  half  of  the  cycle,  the  potential  across  the  diode  becomes  negative, 
and  some  of  the  injected  charge  initially  stored  in  the  neutral  region  of  the 
diode  returns  to  the  driving  circuit.  However,  for  devices  where  the  minority 
carrier  lifetime  is  much  longer  that  the  rf  period,  it  can  be  shown  that 
during  the  negative  half  cycle,  the  charge  returning  to  the  external  circuit 
is  much  less  than  that  injected,  and  it  would  take  a  period  corresponding  to 
several  rf  cycles  for  most  of  the  charge  to  leave  the  diode.  Consequently,  on 
the  next  cycle,  when  the  applied  potential  becomes  positive,  the  voltage 
across  the  diode  does  not  reach  as  large  a  positive  potential  as  in  the  first 
cycle.  It  therefore  follows  that  after  a  few  cycles,  when  the  applied 
potential  becomes  positive,  only  a  very  small  amount  of  excess  charge  is 
injected  into  the  diode  on  each  positive  half  cycle.  Consequently,  the 
potential  across  the  diode  can  be  estimated  easily  by  regarding  it  as  a  simple 
capacitor  and  neglecting  the  conduction  current.  It  is  therefore  convenient 
to  model  the  diode  as  a  simple  capacitor  to  determine  the  rf  voltage  present 
at  the  plate,  and  from  this  information  calculate  the  charge  that  recombines 
in  the  diode. 

The  rf  voltage  across  the  diode  is 

V(t,z)  =  y  V  cos  wt  +  Y  V  cos  (wt  -  kz)  (1) 

W  W  u  o 

where 
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(2) 


Cft  o  o  o  -1/2 

\  -  Hi ♦  r-  -  “2lc„)2  ♦  (~Rpc0)23 

OX  H 

and 

Ya  ■  t(l  ♦  ^)2  *  (“RaC0)2f1/2  (3) 

3 

where 


C 


o 


Cd  +  Cd 


(4) 


Only  the  magnitudes  of  yw  and  Yfl  have  been  used  here  since  any  phase 

change  associated  with  the  phases  of  y  and  y,  can  be  removed  through 

appropriate  definitions  of  the  zeros  of  t  and  z  .  Using  the  values  in 

Table  I,  y  =  .14  .  Equation  (1)  can  be  rewritten  as 
w 

V(t,z)  =  V  cos  (wt  -  <J>)  (5) 


where 


=  C(TwVw  +  TaVa  cos  kz)l 


?  1/2 

(rava  sin  kz)  j 


and 


♦ 


\Vw  +  Va  C0S  kz/ 


(6) 


(7) 
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To  first  order  in  Y  V  h  V  ,  Eqs.  (6)  and  (7)  become 

a  d  W  W 

V  =  Vw  +  TaVa  cos  kz  (8) 


and 


„  Va 

Y  V 
w  w 


sin  kz 


(9) 


These  expressions  will  be  used  to  calculate  the  recombined  charge.  Note 

that  we  expect  the  theory  to  break  down  for  small  plate  voltages  where 

Equations  (8)  and  (9)  are  not  valid  and  also  for  small  times  (i.e.,  the  first 

few  rf  cycles)  where  the  calculation  for  Y  and  y^  is  not  valid  since  a 

w  a 

significant  portion  of  the  diode  current  is  charge  flowing  into  the  neutral 
region. 

The  minority  carrier  density  pn(x,t)  in  the  diode  obeys  the  diffusion 
equation: 


(10) 


where  x  is  the  distance  into  the  neutral  n  region,  Dp  is  the  diffusion 

coefficient  for  the  holes,  and  t  is  the  recombination  time  of  the  holes. 

P 

It  can  be  shown  that  the  total  excess  minority  carrier  charge  in  the  neutral 
region  QD  is  given  by  the  solution  of  Eq.  (9)  in  the  form 
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(11) 


P„(T) 

(t  :,)*«  " 

where  Pn(T)  is  the  excess  minority  carrier  density  at  a  time  t  at  the 
plane  x  =  0  .  The  charge  that  recombines  in  the  diode  (Qr)  is  given  by  the 
differential  equation 


dQ^  _  Qjj 
dt  'p 


Oust  as  in  a  sample  and  hold  circuit,  the  charge  which  recombines  in  the 

diode  Qd  is  charge  stored  on  the  external  capacitor.  This  produces  a  change 
K 

in  voltage  across  the  diode 


VR  *  -V(C0 


+  c  ) 

ox7 


(13) 


Since 


p„(t>  ■ 


evn(t) 


rno 


(14) 


where  $  ~  q/(kT)  and  pn0  is  the  equilibrium  hole  concentration  on  the 
n  side  of  the  junction,  it  follows  that 


3Vsi  n  (tuT+«t>) 
e 


•<v«w 

e 


~(t-T)/T 


(t  -  t) 


p 

T/7 


dr 

(15) 


We  shall  consider  the  case  where  the  minority  carrier  lifetime  x^  is  much 
longer  than  an  rf  period.  Then,  the  most  rapid  variation  occurs  in  the 
exponential  term  containing  the  rf  dependence,  and  we  obtain  the  approximate 


result 


Q  “OP 
sp  HHno 


-fV(Cox+Co>  ~(t  )/Tp 

[I  (ev)e  R  ox  0  -1]  £ - m  dx  (16) 

°  .  (t  -  x)1/z. 


where  I0  is  the  modified  Bessel  function  of  order  zero.  For  slow  charging, 
the  first  term  in  brackets  is  a  slowly  varying  function,  and  the  second  term 
in  brackets  is  only  significant  for  x  «  t  .  Thus,  evaluating  the  first 
term  at  t  =  x 


Q  “  qp 
xp  MKno 


€ 


CI0(6V)e 


-SQR(t)/(C0x4C0) 


-»/ 


t  -<l-T>/Tp 
re  H 


(t  -  x)J 


Qp  “  ^no  VVp  Cl0(eV)e~P9R/{C0X+C0)  “  1]  erf  (t/xp)1/2  (18) 


where 


2  ?  -u2 

erf  x  =  —  J  e  du 


Using  Eq.  (11),  we  obtain  a  simple  differential  equation  for  Qn 


<K).  [oT  -eq„/(c  +c  )  /t  \ 

-  qpno  I0(6V)  (e  .l)erf|_j 


This  has  the  solution 


in 


I.  -  e 


-e 


«lR/(c0X«0) 


6qp-n-0  ■  J5~Tf(t) 

2<?ox+Co>  V  P  P 


where 


■  c-iHef 

r»  n  *  rv 


and  q  is  the  charge  previously  stored  in  the  array.  In  some  of  the 
KO 

figures  in  the  next  section,  the  writing  signal  is  on  for  a  longer  period  than 
the  acoustic  signal.  Consequently,  we  need  to  calculate  the  charge  stored 
when  only  the  writing  signal  is  present  QRg  and  then  calculate  the  final 
charge  QR  after  both  the  writing  and  acoustic  signals  are  present.  Thus, 

q  =  1°* ..L-P-  in  [i  (bv)  -  e"Af(t)  (I  (3V)  -  eeQRo/(Cox+Co^-j  (22) 


where 
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A  * 


(23) 


2<Cox+Co> 

When  QrQ  *  0  ,  Eq.  (22)  is  equivalent  to  the  result  obtained  in  ref.  6  for 
the  air-gap  correlator. 

The  derivation  given  here  has  neglected  the  recombination  of  the  injected 
charge  ,  which  occurs  after  the  pulse  has  been  turned  off.  This  effect  is 
negligible  for  writing  times  long  compared  to  the  minority  carrier  lifetime 

TP  * 

It  is  interesting  to  note  that  if  a  Schottky  characteristic  had  been 
assumed 


qV0/(kT) 


(24) 


then  the  same  result  as  Eq.  (22)  would  have  been  obtained,  except  that  the 
argument  of  the  exponential  term  would  have  been 


instead  of 


ox 


+  C, 


Ist 


Af(t) 


(25) 


for  t/^  »  1 


Assuming  a  typical  value  for  Schottky  diodes 
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(Is  »  10 amps/cm2)  ,  the  charging  rate  is  15  times  faster  than  for  p-n 
diodes.  However,  the  assumptions  of  the  theory  are  questionable  for  such  fast 
charging  times. 

The  recombined  charge  given  in  Eq.  (22)  is  essentially  the  charge  across 
the  plate  capacitor  when  the  diode  voltage  is  at  its  rf  peak.  However, 
after  the  plate  signal  is  removed,  the  charge  redistributes  between  the  diode 
and  plate  capacitors.  If  the  variable  capacity  of  the  diode  under  reverse 
bias  is  taken  into  account,  then  the  relation  between  the  total  recombined 
charge  Qr  and  the  final  charge  on  the  diode  itself  Qrd  after 
redistribution  is 


Q„(l  ♦ 

U0X 


2C0COX  M1/2  , 

(c  ten  / 

1  ox  o'  yo  ' 


(26) 


where  Q0  is  the  depletion  layer  charge  at  thermal  equilibrium.  In  our  case, 
Qr  «  Qq  for  all  of  the  charging  voltages  considered  in  the  next  section. 

In  this  case,  Eq.  (26)  becomes 


QRD 


(27) 


This  is  exactly  what  one  would  calculate  if  the  diode  capacity  were  assumed  to 
be  constant. 

If  we  now  include  the  spatial  variation  of  the  acoustic  wave,  Eq.  (8), 
then  the  redistributed  charge  is 
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Qrd  .  fa  C,  -  e-Af(t)  (I  .  e8<W(Cox+Co));| 

3 


where 


1  =  IoCB(Vw  +  YaVa  cos  kz)]  (29) 

The  charging  time  of  the  diode  in  the  theory  was  taken  to  be  t  .  Since 
the  plate  voltage  is  generally  much  larger  than  the  internal  acoustic  voltage, 
this  charging  time  is  determined  by  the  duration  of  the  write-in  signal  to  the 
plate  tw  .  If  the  acoustic  signal  is  present  at  the  diode  during  the  entire 
charging  time,  then  the  stored  charge  is  given  by  Eq.  (23)  with  t  =  tw  . 
However,  if  the  plate  signal  is  present,  either  before  or  after  the  acoustic 
signal  is  present,  then  there  is  an  additional  component  of  charge  with  no 
spatial  variation.  Using  the  fact  that  Af (t )  <<  1  ,  and  using  Eq.  (22)  to 
obtain  QR  and  QRo  ,  then  the  following  result  is  obtained 

Q 

Qrd  ■  7  F  +  Af(‘w  -  ‘a)  >0  «Vw>  +  Af<V  'o  C»t\VH  +  V.»| 

(30) 


for  tw  >  t ^  ,  and 


{l  *  Af(tA)  i„  Cc(ywvw  *  rava)3{ 


for  tw  <  t^  . 


(31) 
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This  Is  the  expression  for  the  total  stored  charge.  For  the  read-out 
process,  we  are  only  concerned  with  the  component  of  stored  charge  at  the 
acoustic  spatial  wave  number  k  . 


/* 

Qk(t)  =  1/2 ir J  QRD(t,z)  cos  (kz)  d(kz) 
0 


(32) 


This  integral  is  calculated  numerically  for  the  figures  in  the  next  section. 
However,  a  good  approximate  relation  may  be  found  by  assuming  that  the  second 
and  higher  harmonics  are  all  negligible,  and  writing 


For  tw  >  ty\  , 


Qk(t) 


Qmax  (t)  _pmin(t) 
2 


(33) 


<>k 


5>  ln  1  4  *  -  V  'o  t  V,>  *  «V.  *  V,» 

26  1  ♦  A  f(t„  -  tA)  I0(6YKVW)  ♦  A  f(tA)  I0  [6(rwVw  -  T,V,)]  ^ 


For  tw  <  ta  » 


q  .  Eo1n1.l-Af(tw)1°Ce(Vw.tVan  (35) 

26  1  +  Af(tw)  I0C6(ykVw  -  „aVa)] 

Several  results  are  clear  from  this  expression.  Suppose  the  plate  signal 
is  on  for  a  fixed  time  (tw  a  10  ms  for  example)  and  the  acoustic  signal 
length  t^  is  varied  from  0-10  ms  .  In  Eq.  (34),  the  third  term  in  the 
numerator  and  denominator  is  small  compared  to  the  first  two  terms.  Putting 
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f(t)  -  t/Tp  ,  and  carrying  out  a  Taylor  expansion,  we  see  that  is 
linear  with  respect  to  t^  ,  except  in  the  region  where  t^  **  tw  .  This  is 
basically  because  the  presence  of  fixed  charge  due  to  the  applied  plate 
potential  over  the  period  of  10  ms  does  not  affect  the  spatial  modulation 
due  to  the  acoustic  term. 

Suppose  instead  that  the  duration  of  the  acoustic  signal  is  fixed  (10  ms 
for  example),  and  the  duration  of  the  signal  on  the  plate  is  varied 
from  0  -  10  ms  .  Using  Eq.  (35),  it  can  be  seen  that  the  second  term  in  the 
numerator  and  denominator  is  initially  smaller  than  1  ,  and  later  much  larger 
than  1  .  The  net  effect  is  a  very  nonlinear  characteristic.  Put  another 
way,  the  variation  of  the  plate  signal  also  affects  the  spatial  modulation  of 
the  stored  charge  due  to  the  acoustic  signal,  as  the  uniform  charge  stored  is 
no  longer  invariant  with  time. 

From  the  form  of  Eqs.  (30)  and  (31),  (or  even  from  the  approximate  result 
of  Eqs.  (34)  and  (35),  it  is  not  obvious  that  the  output  should  be  linear  with 
respect  to  Vfl  >  except  for  the  case  when  3ywVw  »  1  and  the  expressions 
for  Iq  and  the  logarithm  may  be  expanded  and  a  linear  dependence  on  Va  is 
obtained.  Using  computer  calculations,  it  will  be  shown  in  the  next  section 
that  the  theoretical  dependence  of  Qk  on  Va  is  linear,  and  consequently, 
we  shall  express  our  relations  for  Qk  =  h(Vw)Va  .  For  large  Vw  , 

Eqs.  (34)  and  (35)  are  independent  of  Vw  . 

(36) 

<!„  2VM 

where  Vfoj  is  the  built-in  voltage  of  the  diode. 
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The  Readout 


Plate  to  acoustic  read-out  is  treated  in  this  section.  In  this  case,  an 
rf  modulated  signal  is  applied  to  the  plate  port  and  interacts  with  the 
spatially-modulated  capacity  of  the  diode  array  to  generate  two  acoustic  waves 
which  travel  towards  each  transducer.  We  shall  assume  that  the  amplitude  of 
the  stored  charge  is  constant  across  the  diode  array  and  that  the  read  signal 
is  a  square  pulse  of  duration  t^  .  For  more  complicated  signals,  it  is 
straightforward  to  show  that  the  envelopes  of  the  two  acoustic  waves  are  the 
correlation  and  convolution  of  the  stored  signal  and  the  read  signal. 

We  shall  assume  that  the  read  voltage  is  never  large  enough  to  forward 
bias  the  diodes;  if  the  read  voltage  were  that  large,  the  stored  charge 
pattern  would  be  altered  during  the  read-out  process,  and  the  analysis  is  more 
complicated. 

The  voltage  at  the  diode  due  to  the  read  voltage  will  be  calculated 
first.  Then  a  normal  mode  approach  will  be  used  to  calculate  the  amplitude  of 
the  excited  wave. 

Using  the  circuit  model  in  Fig.  2c,  the  relation  between  the  plate 
voltage  and  the  diode  voltage  is 


V„(z> 


1  +  juRC0(z) 

1  +  jwRC0(z)  +  C0(Z)/C0X 


where  CQ(z)  =  (^(z)  +  .  Oue  to  the  plate  lead  inductance,  the  voltage 

at  the  plate  Vp  is  different  from  the  externally  measured  read  voltage 


(38) 


iy  =  °ivr!  • 


and 


a  *= 


1  +  a>2R2C2 


(1  - 


.2LCp)2 


1/2 


(39) 


where  Cp  is  the  series  combination  of  Cox  and  C0  . 

The  diode  capacity,  (^(z)  is  easily  found  in  terms  of  the  stored  charge 
Qr(z)  • 


1 


Qr(z) 


C  +  C  Q 
ox  o  yo 


(40) 


This  expression  is  valid  to  first  order  in  Qr/Q0  .  Vlhen  Eq.  (40)  is 
substituted  into  Eq.  (37),  the  spatially-varying  component  of  Vd(z)  is  found 
to  be 


'ox 


C  +  C 
o  ox 


1  +  if  +  >RCox 


(41) 


The  amplitude  of  the  surface  acoustic  wave  which  is  excited  by  this 
potential  can  be  found  using  the  normal  mode  analysis  oF  Auld  and  Kino.*5 
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where  k  Is  the  wave  number  of  the  acoustic  wave,  L  is  the  total  length  of 
the  diode  array,  6  is  the  diode  spacing,  l  is  the  diode  period,  and  t^/T 
is  the  ratio  of  the  durations  of  the  read  and  stored  signals.  For  t^/T  >  1  , 
the  ratio  should  be  replaced  by  1. 

Correlator  Output 

To  express  the  output  in  terms  of  external  parameters,  we  need  to  relate 
the  internal  acoustic  voltages  Va  and  $a  to  the  external  input  and  output 
electrical  voltages  V/^  and  Vout  ,  respectively. 


where  IL  is  the  insertion  loss  in  dB  ,  Za  is  the  acoustic  impedance,  and 
Rs  is  the  source  impedance. 

The  external  output  voltage  of  the  correlator  can  now  be  calculated  in 
terms  of  external  parameters 

¥«,t  ■  »<»„>  &  Va  <«> 

where 
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IV.  Results  and  Discussion 
Correlator  Efficiency 

The  dependence  of  the  output  on  and  Vr  is  linear;  however,  the 
output  saturates  as  Vw  is  increased  (Eq.  35).  Consequently,  the  terminal 
efficiency  Ft  shall  be  defined  as  the  output  power  divided  by  the  available 
powers  associated  with  the  two  linear  parameters.  The  terminal  efficiency 
(dBm)  is 


(47) 

(48) 

where 

hmax  =  (49) 

bi 

and  is  the  power  loss  which  results  when  the  plate  is  not  tuned 

to  50  ft  .  Assuming  the  amplifier  is  well  represented  by  an  ideal  voltage 

source  in  series  with  a  source  resistor  Rc 

5  > 


10  log 


out  ' 
PA  PRE 


Ft  =  20  log  (.325  hmax  yr) 
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(50) 


r  (i  -  App2  +  Ay 

[u  -  Acp)2  +  w2C2(R  +  Rs)2 


Using  the  parameters  given  in  Table  1,  y^  =  .5  .  Experimentally,  it  is  found 
that  Yjj  is  closer  to  1  ,  although  the  exact  value  depends  on  the  type  of 
amplifier  and  power  levels  used.  The  maximum  terminal  efficiency  is  therefore 


,Rco*cox 

r<i  -  Ac/  +  «prp)2]1/z 


We  shall  examine  each  of  these  terms,  and  evaluate  the  terms  using  the 
values  listed  in  Table  I. 

The  first  term  is  essentially  the  diode  built  in  voltage  normalized  to 
dBm  .  A  50  n  load  on  the  acoustic  port  was  assumed.  This  term  has  the 
value  .04  .  The  second  term  is  the  reduction  in  output  due  to  the  use  of 
discrete  diodes.  In  our  case,  the  diode  width  and  spacing  are  both  4  u  ,  so 
this  term  has  the  value  .24  .  Decreasing  the  diode  widths  and  spacing  while 
keeping  the  width  equal  to  the  spacing  only  negligibly  increases  this  term  to 
.25  .  However,  if  the  diode  spacing  is  reduced,  this  term  can  be  increased 
to  .5  .  The  limit  of  .5  is  due  to  the  fact  that  either  the  forward  or 
backward  travelling  wave  is  used,  not  both. 
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The  third  term  is  due  to  the  insertion  loss  and  has  the  value  .16  .  The 
fourth  term  is  a  modified  coupling  coefficient,  and  has  the  value  5.83  .  The 


last  term  results  from  the  conversion  of  a  charge  variation  to  a  capacitance 
variation,  and  the  conversion  of  plate  voltage  Vp  to  Vrjt  .  This  term  has 
the  value  .069  .  Decreasing  the  length  of  the  plate  wire  so  that  the  proper 
tuning  is  obtained  only  slightly  increases  this  term  to  .071  . 

If  an  epitaxial  layer  on  a  low  resistivity  substrate  is  used  to  decrease 
R  ,  then  the  first  term  increases  to  .12  and  increases  the  fourth  term  to 
.077  . 

The  first  and  last  terms  are  monotonically  increasing  functions  of  CQ 
and  1/R  and  consequently  a  low  resistivity  substrate  is  suggested.  However, 
since  the  insertion  loss  also  increases  with  decreasing  resistivity,  there  is 
an  optimum  value.  This  trade-off  is  shown  in  Fig.  4.  For  our  device 
configuration,  the  optimum  resistivity  is  5  x  10^  cm~3  . 

Using  the  values  given  above,  we  predict  a  correlator  efficiency  of 
-64  dBm  .  We  experimentally  observed  an  efficiency  of  -66  dBm  .  Note  that 
if  all  of  the  above  suggestions  are  incorporated,  the  efficiency  can  in 
principle  be  improved  by  17  dB  . 

V.  Experimental  Results 

A  series  of  experimental  results  were  compared  with  theory  using  the 
device  described  in  the  introduction.  Most  of  the  experimental  measurements 
were  carried  out  using  Hewlett  Packard  calibrated  attenuators  and  a  Hewlett 
Packard  oscilloscope.  The  accuracy  of  such  measurements  is,  at  best,  1  dB  , 
and  is  subject  to  further  errors  at  low  signal  levels  due  to  noise. 


In  each  of  the  figures,  the  theoretical  results,  Eqs.  (30)  to  (32),  are 
represented  by  a  solid  line.  The  approximate  results,  Eqs.  (34)  and  (35),  are 
represented  by  a  dashed  line.  The  same  parameters  were  used  in  the 
calculations  for  all  the  figures. 


Time  Dependence 

The  dependence  of  the  output  on  the  acoustic  signal  duration  when  the 
plate  signal  duration  is  kept  constant  at  time  t  =  10  ms  is  shown  in 
Fig.  5.  As  anticipated  in  Section  III,  the  dependence  is  linear  except 
when  t^  **  t  .  In  this  case,  both  the  experimental  and  theoretical  results 
indicate  a  slight  deviation  from  a  linear  curve.  It  can  be  seen  that  the 
approximate  formula  is  accurate  enough  for  most  purposes. 

The  dependence  of  the  output  on  the  write  signal  duration  for  short 
times  (<  1  us)  is  shown  in  Fig.  6.  As  expected,  the  agreement  is  good 
except  for  the  initial  charging  period  when  our  expressions  for  y^  and  y^ 
are  not  valid. 

For  the  slow  charging  input  correlation  mode  (Fig.  7),  the  charging 
dependence  is  similar  to  the  dependence  for  fast  charging,  except  that  the 
time  scale  is  larger  by  a  factor  of  roughly  50,000  .  The  effect  of  having 
the  write  signal  on  for  times  longer  than  the  acoustic  signal  (4  ms)  is  also 
shown  in  Fig.  7.  The  output  decreases  due  to  the  increase  in  the  component  of 
siorrd  charge  with  no  spatial  variation. 

The  dependence  of  the  output  on  the  read  signal  duration  was 
'«t  i ]  1  y  and  theoretically  found  to  be  linear. 


The  dependence  of  the  correlator  output  on  the  number  of  times  a  given 
correlation  is  written  into  the  diode  array  is  shown  in  Fig.  8.  The  form  of 
the  dependence  is  important  for  applications  such  as  adaptive  filtering^ 
where  multiple  correlations  are  added  together.  For  most  applications,  it 
would  be  desirable  to  have  a  linear  dependence,  however,  as  shown  in  Fig.  8, 
this  is  not  the  case.  For  adaptive  filter  applications,^  the  feedback  gain 
should  be  decreased  as  a  function  of  time  (equivalently,  the  number  of  writing 
pulses)  and  consequently  the  dependence  shown  in  Fig.  8  is  desirable. 

Amplitude  Dependence 

The  dependence  of  the  output  amplitude  on  the  acoustic  input  signal 
amplitude  is  shown  in  Fig.  9.  It  can  be  seen  that  the  output  is  very  nearly 
linear  over  a  range  of  35  dB  . 

The  dependence  of  the  output  on  the  writing  signal  on  the  plate  is  shown 

in  Fig.  10  for  four  different  read  voltages.  It  can  be  seen  that  the 

agreement  between  experimental  and  theoretical  results  is  good  except  for 

small  write  voltages  and  large  read  voltages  where  our  approximations  are  no 

longer  valid.  The  charging  period  is  very  long  (t  =  t  =10  ms)  in 

P  a 

Fig.  10  and  the  output  saturates  for  voltages  greater  than  2 V  .  For  short 
charging  periods  (.2  us  ),  as  shown  in  Fig.  11,  the  output  does  not  saturate 
until  the  write  signal  reaches  3V  .  The  error  in  the  write  amplitude  range 
between  2  and  3V  is  due  to  the  fact  that  the  I -V  characteristic  of  the 
diode  follows  an  exponential  curve  of  the  form  exp  qV/kT  .  Hence,  very 
slight  errors  in  estimating  V  show  up  as  large  errors  in  the  curve. 
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The  dependence  of  output  voltage  on  the  read-out  voltage  is  shown  in 
Fig.  12  for  three  different  acoustic  voltages.  It  can  be  seen  that  the 
agreement  with  theory  is  good  except  for  large  read-out  voltages  where  the 
variation  in  diode  capacitance  is  no  longer  linear  with  voltage  as 
approximated  in  the  theory. 

VI.  Conclusions 

We  have  given  here  a  comprehensive  quantitative  theory  which  predicts  the 
output  of  a  monolithic  acoustic  surface  wave  correlator  as  a  function  of  the 
important  external  variables.  An  extensive  set  of  experimental  measurements 
confirms  the  predictions  of  the  theory.  The  theory  can  therefore  be  employed 
for  accurate  design  of  monolithic  acoustic  surface  wave  storage  correlators. 
Suggestions  are  given  for  improving  the  efficiency  of  the  device  by  as  much 
as  17  dB  . 
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FIGURE  CAPTIONS 


1.  Schematic  drawing  of  the  monolithic  storage  correlator.  . 

2.  Circuit  models:  (a)  writing  process;  (b)  diode  model;  and  (c)  reading 
process. 

3.  Experimental  measurement  {•)  of  the  top  plate  impedance  as  a  function  of 
frequency  and  theoretical  fit  { -)  to  the  data  using  the  model  of 

Fig.  2a. 

4.  The  effect  of  doping  level  on  the  acoustic  wave  attenuation  correlation 
efficiency  and  total  correlation  efficiency,  including  attenuation  for 
short  pulses. 

5.  Dependence  of  correlator  output  on  the  duration  of  the  acoustic 
signal. 

6.  Dependence  of  correlator  output  on  the  duration  of  the  write  signal  for 
short  pulses  (0-1  ns). 

7.  Dependence  of  correlator  output  on  the  duration  of  the  write  signal  for 
long  pulses  (0-10  ms). 

8.  Dependence  of  correlator  output  on  the  number  of  correlations. 

9.  Dependence  of  correlator  output  on  acoustic  power. 

10.  Dependence  of  correlator  output  on  writing  signal  voltage  for  four 


reading  voltages  and  10  ms  input  correlation  times. 

11;  Dependence  of  the  correlator  output  on  writing  signal  voltage  for  two 
acoustic  power  levels  and  .2  ps  input  correlation  times. 


12.  Dependence  of  output  correlation  on  the  readout  voltage  for  three 
different  acoustic  voltages. 


TABLE  I:  Parameters  of  the  Storage  Correlator  Used  jn  the  Experiments 

Measured  Parameters 


L 

Length  of  diode  array 

1.1  cm 

1 

Diode  periodicity 

8.0  iim 

6 

Diode  spacing 

4.0  pm 

X 

Acoustic  wavelength 

Suhstrate  thickness 

Si02  thickness 

32.0  pm 

350.0  pm 
.1  pm 

ZnO  thickness 

1.6  pm 

Av/v 

Coupling  coefficient 

.0027  (.0036*) 

"p 

Substrate  resistance 

5.0  ft  (5.2*) 

i 

Top  plate  lead  inductance 

22  nH  (20  nH*) 

cp 

Series  capacitance  of  ZnO  and  Si02 

Top  plate  area 

Type  of  substrate 

Calculated  Parameters 

364.0  pF  (326  pF*) 
.11  cnr 

N  -type 

u 

Angular  frequency 

7.8  x  103  rad/s 

Cox 

ZnO  and  Si02  capacity 

5850  pF/cm^ 

^a 

Acoustic  capacity 

7950  pF/crn2 

cd 

Diode  capacity 

6030  pF/cm^ 

Co 

Silicon  capacity  (C^  +  C^') 

6030  pF/cm3 

Np 

Substrate  doping  density 

Estimated  Parameters 

3.5  x  lO1'1  cm'3 

T 

P 

Minority  carrier  lifetime 

.1  PS 

Cd 

Capacity  between  the  diodes 

6030  pF/cm2 

*  Theoretical  calculation 


RELATIVE  CORRELATION  EFFICIENCY  (dB) 


DOPING  DENSITY  (cm 


P0  =  5  dBm 


ACOUSTIC  SIGNAL  DURATION  (ms) 


CORRELATION  OUTPUT  (mV) 


WRITING  SIGNAL  DURATION  (ycis) 


CORRELATION  OUTPUT  (ARBITRARY  UNITS) 


CORRELATION  OUTPUT  VOLTAGE  (mV) 


3 


5 


APPENDIX  C 


DISPERSION  EFFECTS 


1.  Introduction 

Any  SAW  device  that  utilizes  a  layered  substrate  -will  be  dispersive  i.  e. 
the  velocity  of  propagation  will  depend  on  the  frequency  of  the  wave.  The 
SAW  velocity  of  the  material  used  for  the  substrate  will  in  general  be  different 
than  the  SAW  velocity  of  the  material  used  for  the  layer.  Consequently, 
the  SAW  velocity  of  the  layer/substrate  system  will  depend  on  the  relative 
amounts  of  energy  in  the  layer  and  in  the  substrate.  As  the  frequency  of  the 
wave  decreases,  the  decay  depth  of  the  surface  wave  increases  and  the  SAW 
velocity  approaches  the  SAW  velocity  of  the  substrate.  Dispersion  will  thus 
be  minimized  when  the  SAW  velocities  of  the  substrate  and  layer  are  almost 
equal.  Thus  systems  like  ZnO  (vr  =  2639)  on  GaAs  (y«  =  2600)  are  not 
very  dispersive  while  systems  like  ZnO  on  silicon  (vjj  =  4920)  or  ZnO  on 
sapphire  (vr  =  5600)  are  very  dispersive. 

Dispersion  is  of  course  not  a  problem  in  narrow  band  systems  since  the 
velocity  does  not  change  very  much  over  the  frequency  range.  However,  dis¬ 
persion  causes  phase  distortions  of  the  wave,  and  in  broadband  systems,  dis¬ 
persion  can  limit  the  usable  bandwidth  to  much  less  than  the  IDT  bandwidth. 


In  this  section  we  shall  present  results  on  the  variation  of  group  velocity 
with  frequency  and  compare  the  amount  of  dispersion  in  the  different  systems 
discussed  in  the  previous  section.  Then  we  shall  analytically  examine  the 
effects  of  dispersion  on  correlation  and  convolution. 

2.  Group  Velocity 

The  dependences  of  the  group  and  phase  velocities  on  normalized  wave 
number  hk  for  the  Rayleigh  mode  in  ZnO  on  (100)Si  and  (lll)Si  are  given  in 
Fig.  2.1.  The  position  of  the  two  peaks  in  Av/v  are  indicated  in  the  figure. 
Note  that  the  second  Av/v  peak  occurs  at  the  optimum  value  of  hk,  namely 
the  point  where  the  group  velocity  goes  through  a  minimum.  The  curves  for 
(100)Si  and  (lll)Si  have  the  same  amount  of  dispersion  at  the  position  of  the 
second  Av/v  peak;  however,  at  the  first  Av/v  peak  the  slope  of  the  group 
velocity  curve  for  the  (lll)Si  system  is  half  the  slope  of  the  (100)Si  system. 

The  same  data  for  Sezawa  waves  on  (100)Si  is  presented  in  Fig.  2.2.  In 
this  case,  the  horizontal  axis  is  multiplied  by  the  phase  velocity  to  obtain 
hu  which  unlike  hk  is  an  experimentally  measurable  quantity.  This  is  done 
so  that  the  experimental  data  could  be  included.  It  can  be  seen  that  the 
agreement  is  very  good. 

These  systems  can  be  numerically  compared  by  tabulating  the  slope 
of  the  group  velocity  curve  at  the  desired  operating  point  (Table  1).  This 
parameter  will  be  used  for  the  analysis  in  the  next  section;  however,  a  better 
physical  grasp  of  the  situation  can  be  obtained  by  calculating  how  much 
spreading  occurs  to  a  .04  pa  pulse  after  it  has  traveled  1  cm  (column  2,  Table 
1)  or  after  it  has  traveled  1  pa  (column  3,  Table  1).  The  magnitude  of  the 
dispersion  problem  is  now  apparent.  For  instance,  suppose  it  was  desired 
to  delay  a  digital  code  (25  MHz  bit  rate)  by  1  pa.  The  output  would  be 
hopelessly  distorted  by  any  of  the  monolithic  systems  considered  here.  The 


maximum  bit  rate  that  could  be  used  with  a  Sezawa  wave  1  pts  delay  line  is 
14  MHz.  Thus,  the  usable  bandwidth  is  entirely  limited  by  dispersion,  not 
the  IDT  bandwidth. 

For  systems  like  ZnO/Si  Sezawa  wave  or  Rayleigh  wave  (first  Av/v 
peak),  the  variation  of  group  velocity  is  linear.  Thus,  a  chirp  transducer  with 
linear  phase  dispersion  can  be  used  to  compensate  for  the  material  dispersion. 
This  will  work  for  any  delay  line  with  a  fixed  delay.  The  technique  is  only 
approximate  for  devices  like  convolvers  and  correlators  where  the  delay  to 
different  points  in  the  interaction  region  varies  from  typically  0  to  4  ps. 
However,  as  we  shall  see  in  the  next  section,  compensation  can  be  obtained 
for  many  systems  utilizing  correlators  or  convolvers. 

3.  Effect  of  Dispersion  on  Correlation  and  Convolution 

It  is  useful  to  qualitatively  consider  the  effects  of  dispersion  on  1)  correlar 
tion  and  input  correlation  in  the  storage  correlator,  and  2)  convolution  in  the 
SAW  convolver.  We  shall  then  analytically  examine  the  effects  of  dispersion 
on  these  processes. 

.  Qualitative  Dispersion  Effects 

We  shall  use  a  chirp  signal  for  our  discussion  since  its  use  has  important 
signal  processing  applications  and  since  the  different  frequency  components 
are  clearly  distinguishable. 

The  correlation  of  a  chirp  is  shown  in  Fig.  2.3.  For  the  ZnO/Si  system, 
the  group  velocity  decreases  with  increasing  frequency  (for  hu  <  104)  so  for 
a  down  chirp,  the  initial  high  frequency  components  travel  the  slowest,  and 
the  chirp  contracts.  This  signal  is  stored  by  applying  a  short  signal  to  the 
topplate.  A  second  signal  is  sent  to  the  same  acoustic  port;  it  undergoes 
the  same  distortion  and  the  correlation  is  obtained  at  the  topplate  port. 


The  output  is  compressed  from  the  initial  width  by  the  full  time  bandwidth 
product  of  the  device.  However,  for  times  before  and  after  the  time  of  chirp 
coincidence,  the  distortions  on  the  two  signals  are  different,  and  thus,  the 
sidelobe  pattern  of  a  dispersive  device  will  be  different  than  the  sidelobe 
pattern  of  a  nondispersive  device.  Thus,  the  width  of  the  output  pulse  is 
not  increased  but  the  sidelobe  pattern  will  change.  This  change  is  usually 
unimportant.  Note  that  these  conclusions  are  true  for  any  compression 
experiment.  The  code  can  be  a  chirp,  PSK  (phase  shift  keyed),  or  frequency 
hopped,  and  the  maximum  bit  rate  allowed  by  the  transducer  bandwidth  can 
be  used. 

For  input  correlation  in  a  device  without  dispersion,  two  long  chirps 
are  sent  into  the  device,  and  a  narrow  peak  builds  up  in  the  center  of  the 
interaction  region.  Much  larger  compression  ratios  can  be  obtained  using  the 
input  correlation  mode  (30,000  for  a  ZnO/Si  correlator [1])  since  the  duration 
of  the  codes  is  only  limited  by  the  storage  time  of  the  diodes.  However,  the 
dynamic  range  is  limited  because  all  of  the  charge  is  stored  in  just  a  few 
diodes,  and  the  full  interaction  region  is  not  being  utilized. 

Consider  now  a  dispersive  correlator.  Only  one  of  the  two  chirp  signals 
is  distorted,  unlike  the  regular  correlation  mode.  The  effect  of  distorting  only 
one  chirp  is  that  as  the  correlation  proceeds,  the  position  of  the  narrow  pulse 
changes  (it  shifts  from  left  to  right  in  Fig.  2.3).  The  total  peak  shift  is  given 
by  the  change  in  group  velocity  over  the  chirp  frequency  range  times  the  time 
delay  to  the  center  of  the  interaction  region.  This  width  is  .45  cm  for  a  typical 
Sezawa  wave  device.  Fortunately,  the  phase  relations  in  the  broadened  pulse 
are  such  that  when  an  rf  read  pulse  is  applied  to  the  acoustic  port,  the  full 
compression  should  be  obtained.  Furthermore,  the  dynamic  range  should  be 
better  than  for  a  nondispersive  device  since  much  more  of  the  interaction 
region  is  utilized  in  the  dispersive  device.  Thus,  dispersion  does  not  degrade 


the  storage  correlator  performance  in  compression  experiments. 

In  a  convolver,  two  counterpropagating  waves  are  sent  into  the  device 
(Fig.  2.4).  The  effect  of  dispersion  is  to  compress  one  chirp  and  expand  the 
other.  Thus,  the  output  pulse  width  will  be  broadened  by  dispersion.  Note, 
however,  that  in  a  device  with  a  linear  change  in  group  velocity  that  the  full 
compression  can  be  obtained  by  appropriately  changing  one  of  the  slopes, 
either  by  using  a  chirp  transducer  or  by  using  two  chirp  generators  with 
different  slopes. 

Analytic  Effects  of  Dispersion 

Consider  first  a  correlator  operating  in  the  regular  correlation  mode.  A 
signal  ai(t)cosu/ot  is  sent  into  the  acoustic  port.  If  the  Fourier  transform  of 

this  signal  is  Ai(w)  i.  e. 

1 

ai(t)cosw0t= -  /  dwAi(w)eJwt  2.49 

After  this  signal  has  propagated  for  a  time  to,  a  short  pulse  is  applied  to  the 
topplate,  and  the  signal  q(z )  is  stored 

q{z)  =  —  f  2.50 

V2irJ-co 

For  acoustic  to  plate  readout,  a  second  signal  O2(t)cosu/0t  with  Fourier  trans¬ 
form  A2(w)  is  sent  to  the  acoustic  port,  and  the  time  and  spatial  dependence 
of  the  propagating  signal  is 

o2(«,  t)  =  —  f  dw,Aa(w,)e,(w<— 2.51 

V2kj-°° 

The  output  is 

V(t)  =  I  dzaa(z,t)q(z)  2.52 

Jo 


V{t)  =  i-  J  dwdw'A1(w)A2(w')  I 

V(t)  =  J  dwdw,A1(w)A2(a/)^-«<o)tf(ife(a;)  _  *(„*)) 


Since[2] 


/  ?(*)^(/(a:))da:  =  £ 

J— oo  .  irf/i 


Z  l*i 


where  zt-  are  the  zeros  of  f(x),  and  since 


2.53 

2.54 


2.55 


1 

dk 

dw  ui  —  u1 


2.56 


then 

V(t)  =  -^f_Qo  dw|y<?(w)|A1(w)A2(-w)  2.57 

For  no  dispersion,  VGf(w)is  a  constant  (Vgo)  and  the  output  is  just  the  cor¬ 
relation  of  the  two  inputs: 

r 

Vo(t)  =  Vgo  /  dwA1(w)A2(— w)  =  VGfo[ai(^)coso/0t]  ®  [a2(i)cosw0^] 

J — oo 

2.58 

where  (g)  stands  for  correlation.  When  dispersion  is  present,  then  the  output 
is  the  convolution  of  1)  the  correlation  output  in  the  absence  of  dispersion 
Vo(£)  with  2)  the  Fourier  transform  of  the  group  velocity 


V(t)  =  V„*F(|V0(u>)|) 


For  the  case  of  linear  dispersion,  this  transform  is 

F(|Vc(w)|)  =  £vO0(^)s 

T  r 


where 


r 


1  dVg 
V<jo  du> 


2.59 


2.60 


2.61 

v/ 
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Thus,  +he  output  is  broadened;  however,  for  any  of  the  monolithic  devices 
discussed  here,  r  is  so  small  that  the  broadening  is  not  detectable.  For 
instance,  for  a  ZnO/Si  Sezawa  wave  device  at  165  MHz,  r  =  5xl0— 14«. 

For  input  correlation,  the  analysis  is  the  same  except  that  the  stored 
charge  is 

q(z)  =  J  A1(w)A2(w)e-Jw^v'cVw9r*Vo°  2.62 


where  r  is  defined  in  Eq.  2.61.  If  r  is  nonzero,  then  there  is  a  nonzero 
distorting  term.  If  a  delta  function  read  pulse  is  assumed,  then  Eq.  2.58  is 
obtained.  Thus,  even  if  r  is  nonzero,  the  correlation  output  is  not  significantly 
broadened  by  dispersion. 

The  case  of  dispersion  in  a  convolver  has  been  treated  by  Morgan [3]. 
The  analysis  is  similar  to  the  correlation  analysis.  The  two  traveling  waves 
are 


ai{z,t) 


1 

sfa 


/oo 

-OO 


2.63 


a2(z,t)  =  —  f  du'A2{ 


The  output  is 


V(t)  =  f  dzai{z,t)a2(z,t) 

Jo 

/OO 

du\VGo  ( w)  |  Ai  (w)  A2  {u)e^2wt~  k^L 

-inf  tv 


•infty 

Using  the  expansion  for  k 


2.64 


2.65 

2.66 


w  .  d  ,  1  .  2  w  .  r  2 

- - - ( — tar  - - a r 

Vgo  dwVa  VG0  Vbo 


2.67 


and  changing  the  transform  variable  to  t\  =  2t  —  L/Vqo  yields 
V{t)  =  ^  f  dwA1(w)A2(w)cJwtle-jXrw8/vbo 

2  J — oo 


2.68 
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The  important  difference  between  this  expression  and  the  expression  we 
obtained  for  correlation  is  the  rapidly  varying  phase  term  exp( — jlM2r/Vao). 
For  correlation,  the  total  effect  of  dispersion  entered  through  the  frequency 
variation  of  Vq(u).  For  convolution,  this  relatively  slow  variation  can  be 
neglected.  Since  the  Fourier  transform  of  exp( — aw2)  is  1/a  exp( — t2/a)  then 
the  convolution  output  is  convolved  with  the  distorting  factor 

p(c-i(Lw9r/V<so))  Z ieji*/r'2  2.69 

where  ti  =  y/Lr/Vao  .  For  a  Sezawa  wave,  r*  =  .03 pis.  Consequently,  the 
maximum  usable  bandwidth  is  on  the  order  of  30  MHz. 


.  REFERENCES 

[1]  H.  C.  Tuan  and  G.  S.  Kino,  El.  Lett. 

[2]  R.  Bracewell,  "The  Fourier  Transform  and  Its  Applications",  McGraw- 
Hill,  New  York,  95(1965). 

[3]  D.  P.  Morgan,  "Effect  of  Dispersion  in  Surface  Acoustic  Wave  Convolvers", 
IEEE  Trans.  Sonics,  Ultrasonics,  SU-22,  274(1975). 


List  of  Figures 

[1]  Dependence  of  group  and  phase  velocities  on  normalized  wave  number 
hk  for  the  first  Rayleigh  mode  on  (100)Si  and  (lll)Si. 

[2]  Dependence  of  group  and  phase  velocities  on  normalized  frequency  hu 
for  the  Sezawa  mode  (theory  experiment  .). 

[3]  Schematic  drawing  of  chirp  correlation  experiment  a)  regular  correla¬ 
tion  mode,  b)  input  correlation  mode. 

[4]  Schematic  drawing  of  chirp  convolution  (one  chirp  inverted). 

List  of  Tables 

[1]  Dispersion  characteristics  for  several  systems  of  interest. 


9 


APPENDIX  D 


J 
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The  phase  velocities  and  coupling  coefficients  for  Sezawa  waves  in  (001)  ZnO  on  (001)  cut,  (100) 
prop.  Si  and  on  (1 1 1)  cut,  (1 12)  prop.  Si  are  presented.  Experimental  results  from  Sezawa  wave 
delay  lines  are  presented  which  represent  for  the  first  time  bandwidths  in  excess  of  30  MHz  in  a 
thin-film  monolithic  device. 

PACS  numbers:  43.35.Ns,  77.60.  -f  v,  68.60.  +  q,  43.20.Bi 


Acoustic  surface  wave  (ASW)  signal  processing  devices 
such  as  convolvers  and  storage  correlators  can  be  monolith¬ 
ic,  thin-film  devices  (typically  ZnO  on  Si)1  or  hybrid  ar¬ 
rangements  where  a  piezoelectric  substrate  (typically 
LiNbOj)  is  pressed  against  a  semiconductor  (typically  Si).2'3 
The  monolithic  approach  has  the  advantages  of  ruggedness, 
smaller  size,  and  fewer  spurious  signals,  and  the  devices  are 
fabricated  using  standard  planar  processing  and  are  poten¬ 
tially  much  cheaper  to  make.  The  very  important  limitation 
of  monolithic  ASW  devices  has  been  the  inherent  narrow 
fractional  bandwidth  (10-15%)  which  is  a  result  of  the  low 
coupling  coefficient  (du/e<0.004)  attainable  with  thin-film 
Rayleigh  wave  devices.  Larger  first-order  Rayleigh  wave 
coupling  coefficients  (Av/v<0.0l5)  are  attainable  with 
thicker  films  (A— 0.5A.  rather  than  h  —0.07/1 ). 

Armstrong  and  Crampin9  have  presented  theoretical 
calculations  for  the  coupling  coefficient  of  second-order 
Rayleigh  waves,  commonly  called  Sezawa  waves,4  for  (001) 
ZnO  on  (001)  Si.  They  found  a  peak  coupling  coefficient  of 
0.024  at  a  velocity  of 492 1  m/s.  Elliot  and  co-workers7  *  have 
made  Sezawa  wave  signal  processing  devices  with  13% 
bandwidth. 

In  this  letter,  we  present  theoretical  and  experimental 
values  for  the  phase  velocity  and  coupling  coefficient  for 
Sezawa  waves  in  (001)  ZnO  on  (001)  and  (11 1)  Si.  Informa¬ 
tion  is  also  presented  on  the  third-order  Rayleigh  mode  on 
(001 )  Si.  Finally,  results  are  presented  on  Sezawa  wave  delay 
lines  with  fractional  bandwidths  of  19  and  25%  at  a  center 
frequency  of  165  MHz. 

The  dependence  of  phase  velocity  on  the  normalized 
film  thickness  hk  of  the  four  lowest-order  surface  acoustic 
waves  for  a  film  of  (001)  ZnO  on  a  (001)  cut,  (100)  prop.  Si 
substrate  are  shown  in  Fig.  1 .  The  Sezawa  mode  is  cut  off  at 
the  silicon  shear  wave  velocity,  which  is  5800  m/s  for  this 


FIG.  2.  Sezawa  wave  coupling  coefficient  vs  normalized  film  thickness  for 
(001 )  ZnO  on  (1 1 1)  cut,  (1 12)  prop.  Si. 

orientation  and  propagation  direction.  This  mode  is  progres¬ 
sively  more  lossy  as  k  is  decreased  below  the  cutoff  value 
( hkc  =  0.73).  The  cutoff-  value  varies  by  as  much  as  15%, 
depending  on  the  electrical  conditions  at  the  surface  and  the 
film-substrate  interface. 

The  coupling  coefficients  of  the  Sezawa  mode  are 
shown  in  Fig.  2  for  the  four  possible  transducer  configura¬ 
tions.  Note  that  we  predict  coupling  coefficients  as  large  as 
0.028,  which  is  six  times  larger  than  the  first  Rayleigh  wave 
peak  for  this  orientation.4 

The  coupling  coefficient  Av/v  must  be  calculated  at  a 
constant  value  of  to,  not  k  as  in  Armstrong  and  Crampin,9 
since  the  theories4  9  for  transducer  impedance  all  use  Av/v  at 
constant  to.  Our  calculations  for  Av/v  at  constant  k  agree 
with  the  results  of  Armstrong  and  Crampin9  and  yield  a 
value  for  the  peak  coupling  coefficient  which  is  17%  lower 


FIG.  1.  Phase  velocity  vs  nor¬ 
malized  filter  thickness  of  (001) 
ZnO  on  (001)  cut  (100)  prop.  Si. 


FIG.  3.  Phased  velocity  vs  nor¬ 
malized  film  thickness  for  (001) 
ZnO  on  (1 1 1)  cut  (1 12)  prop.  Si. 


FIG.  4.  Sezawa  wave  coupling  coefficient  vs  normalized  film  thickness  for 
(001)  ZnO  on  (1 1 1)  cut.  (1 12)  prop.  Si. 


than  the  correct  value.  Armstrong  and  Crampin  report  a 
phase  velocity  of 492 1  m/s,  which  is  significantly  lower  than 
our  theoretical  and  experimental  results  and  the  experimen¬ 
tal  results  of  Elliot  et  al.1 

The  first-order  Rayleigh  wave  coupling  coefficient  is 
0.0002  at  the  optimum  film  thickness  for  Sezawa  wave  cou¬ 
pling  (hk  =  1.65)  for  transducers  at  the  free  surface  and  a 
metal  shorting  plane  at  the  interface.  Thus  the  spurious  Ray¬ 
leigh  wave  generation  is  50  dB  smaller.  The  coupling  coeffi- 


TABLE I.  Delay-line  parameters. 


Substrate  orientation:  (001) 

L,  s  0.04*  fim 

Number  of  finger  pairs:  4 

L,  =  0.048  fim 

Wavelength:  32  fim 

(,ss0.9gm 

Film  thickness:  8/rm 

C,«36pF 

Beam  width:  1  mm 

C,«40pF 

Pad  capacitance:  0.2  pF 

C,  ss=  9.6  pF 

We  have  built  Sezawa  wave  delay  lines  on  both  the 
(001)  and  (111)  orientations.  The  experimental  agreement 
with  theoretical  values  for  phase  velocity  (Figs.  1  and  3)  and 
coupling  coefficient  (Fig.  2)  is  good. 

Using  the  coupled-resonator  tuning  network  shown  in 
Fig.  5,  a  3-dB  bandwidth  of  31  MHz  with  22.8  dB  insertion 
loss  was  achieved  (Fig.  5).  The  device  parameters  are  given 
in  Table  I.  The  theoretical  plot  of  insertion  loss  was  obtained 
using  the  theoretical  values  for  phase  velocity  and  Au/u  and 
the  measured  values  given  in  Table  I  for  the  device  and  net¬ 
work  parameters.  The  equations  for  transducer  impedance 
and  capacitance  derived  by  Kino  and  Wagers'*  for  balanced 
excitation  were  modified  to  apply  to  the  case  of  unbalanced 
excitation.  Using  a  two-stage  ir  network,10  bandwidths  as 
large  as  42  MHz  were  attained,  although  the  insertion  loss 
was  higher  (35  dB). 

In  conclusion,  theoretical  and  experimental  values  for 
the  phase  velocities  and  coupling  coefficients  are  presented 
here.  Bandwidths  in  excess  of  30  MHz  were  achieved  for  the 
first  time  in  a  monolithic  ZnO-on-Si  device. 


cient  for  Love  waves  to  these  transducers  is  zero,  so  these 
waves  are  not  excited. 

The  maximum  coupling  coefficient  for  third-order 
Rayleigh  waves  was  found  to  be  0.0072  at  hk  —  5.8.  This 
mode  is  cut  off  at  hk  =  3.2  and  is  not  desirable  for  these 
applications  because  of  the  low  coupling  coefficient. 

The  Sezawa  wave  phase  velocity  and  coupling  coeffi¬ 
cient  for  (111)  cut  (112)  prop.  Si  are  shown  in  Figs.  3  and  4, 
respectively.  The  silicon  shear  wave  velocity  is  lower  (5000 
m/s)  for  this  orientation  and  propagation,  and,  consequent¬ 
ly,  the  cutoff  occurs  at  a  higher  value  of  hkc  =  2. 1 5. 


FIG.  S.  Insertion  loss  vs  frequen¬ 
cy  of  a  Sezawa  wave  delay  line  on 
(001)  cut  (100)  prop.  Si. 


The  authors  wish  to  thank  R.  Wagers  and  M.  Drake  for 
their  help  and  guidance  with  the  computer  calculations.  The 
computer  calculations  were  made  using  a  program  devel¬ 
oped  at  Stanford"  and  used  the  crystal  constants  compiled 
by  Auld.12  This  work  was  supported  by  the  Defense  Ad¬ 
vance  Research  Program  Agency  and  monitored  by  the  Of¬ 
fice  of  Naval  Research  under  Contract  No.  N00014-76-C- 
0129. 


'H.C.  Tuan  and  G.S.  Kino.  Appl.  Phys.  Lett.  31,  641  (1977). 

2K.A.  Ingebritsen,  R.A.  Cohn,  and  R.W.  Mountain,  Appl.  Phys.  Lett.  2d. 
596(197$). 

3C.  Maerfeld  and  Ph.  Defrmnould,  IEEE  Ultransonics Symposium  Proceed¬ 
ings — 1975  (IEEE.  New  York,  1975),  p.  20*. 

4G.S.  Kino  and  R.  Wagers.  J.  Appl.  44,  1480(1973). 

’G.A.  Armstrong  and  S.  Crampin,  Electron.  Lett.  9,  322  (1973). 

*K.  Sezawa.  Bull.  Earthquake  Res.  Inst.  3.  I  (1927). 

’J.K.  Elliott,  R.L.  Gunshor,  R.F.  Pierret,  and  AR.  Day,  Appl.  Phys.  Lett. 
32,515(197*). 

"F.C.  Lo,  R.L.  Gunshor,  and  R.F.  Pierret,  Appl.  Phys.  Lett.  34, 725  ( 1979). 
*W.R.  Smith,  H.M.  Gerard,  J.H.  Collins,  T.M.  Reeder,  and  H.J.  Shaw, 
IEEE  Trans.  Microwave  Theory  Tech.  17,  865  (1969). 

"T.M.  Reeder,  W.R.  Shreve,  and  P.L.  Adams,  IEEE  Trans.  Sonics  Ultra¬ 
sonics  19, 466(1972). 

"R.  Wagers,  Ph.D.  thesis,  Stanford  University,  1972  (unpublished). 
I:B.A.  Auld,  Acoustic  Fields  and  Waves  in  Solids  (Wiley-Interscience,  New 
York,  1973),  Vol.  1,  p.  357. 


•07  A**.  Phy*.  LML,  Vol.  30,  No.  10, 15  M«y  1900 


Bowws,  Khuri-Yakub,  and  Kino 


APPENDIX  E 


i  t 
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Abstract 

A  new  analog  adaptive  filter  for  deconvolving  distorted  signals  is 
described  in  this  paper.  The  filter  uses  a  storage  correlator  which 
implements  a  clipped  version  of  the  LMS  algorithm  and  uses  a  special 
iterative  technique  to  achieve  a  fast  convergence.  The  new  filter 
has  a  potential  bandwidth  of  100  MHz  and  would  eventually  handle  pulsed 
signals  of  5  ysec  width.  For  signals  with  time-bandwidth  product  of 
less  than  100,  the  adaptation  time  is  less  than  1  msec  which  allows 
operation  in  real  time  for  most  applications  including  resolution  of 
radar  signals  in  a  cluttered  environment,  removal  of  echoes  from  tele¬ 
vision  signals,  deconvolution  of  distorted  signals  in  nondestructive 
evaluation,  and  also  in  telephony.  The  filter  is  particularly  suited 
for  radar  and  communications  as  it  processes  signals  directly  in  the 
VHF  range. 

Two  experiments  related  to  ghost  suppression  of  a  pulse  and  to  the 
field  of  N'DE  are  described  in  this  paper.  The  results  are  in  good 
agreement  with  computer  simulations  and  show  a  ghost  suppression  of 
15  dB  for  the  first  example  and  a  sidelobe  suppression  of  8  dB  for  a 
transducer  signal.  The  adaptation  time  is  less  than  450  ysec. 
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I.  INTRODUCTION 

Adaptive  filtering*  is  useful  in  removing  distortion  from  signals, 

particularly  when  the  distortion  varies  in  time.  Adaptive  filters  have 

been  used  to  perform  deconvolution  of  a  distorted  echo  pulse  in  an 

2 

acoustic  imaging  system  ,  to  equalize  the  distortion  in  a  telephone 
channel^,  and  to  suppress  an  interfering  signal^. 

Most  adaptive  filters  have  been  implemented  using  digital  techniques. 
The  limitation  of  the  digital  approach  is  the  limited  bandwidth  (typi¬ 
cally  5  MHz)  and  the  practical  limit  on  the  number  of  taps  dictated 
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by  the  complexity  and  power  consumption  . 

An  analog-digital  hybrid  approach  has  been  implemented  using  MOS  LSI 
technology^.  This  has  the  advantage  of  lowering  the  power  consumption 
and  allowing  32  taps  to  be  used  without  undue  external  complexity. 

Large  dynamic  range  was  obtained  with  this  technique  (60  dB) ,  but  the 
bandwidth  was  limited  to  less  than  1  MHz. 

Most  analog  implementations  of  an  adaptive  filter  have  been  made 
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using  CCDs  with  analog  tap  weights  held  in  sample  and  hold  circuits  . 

The  limitations  here  are  the  narrow  bandwidth  and  variation  across  the 
chip  in  gain  and  threshold  levels.  The  alternative  approach  used  at  • 
Hughes  for  implementing  a  wideband  adapting  filter  is  to  employ  tapped 
ASW  filters  with  complex  computer-controlled  systems  for  adjusting  the 
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tap  weignts  . 

We  shall  describe  here  a  relatively  simple  all-analog  approach  to 
adaptive  filtering  which  uses  the  least  mean  squares  (LMS)  algorithm 
to  find  the  optimal  set  of  tap  weights.  An  ASW  monolithic  ZnO/Si  storage 
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correlator  with  8  MHz  bandwidth  and  the  equivalent  of  24  taps  was  used  . 
The  advantages  of  this  approach  are: 
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(1)  fast  iteration  rate  (100  kHz)  which  means  a  short  learning 
time  (~100  ys)  and  good  ability  to  track  time  varying  signals; 

(2)  simple  external  connections; 

(3)  large  potential  bandwidth  (100  MHz); 

(4)  large  potential  number  of  taps  (1000); 

(5)  low  power  consumption  (.1  W); 

(6)  it  is  suited  for  radar  and  communications  systems  because  it 
can  operate  directly  at  the  IF  frequency. 

A  disadvantage  is  that,  at  the  present  time,  the  dynamic  range  is  less 
than  for  digital  systems. 

A  brief  introduction  to  the  LMS  algorithm  is  provided,  followed 
by  a  new  derivation  of  the  LMS  algorithm.  As  the  class  of  applications 
discussed  here  utilizes  continuous,  quasistatic  filters,  a  derivation 
in  the  frequency  domain  rather  than  the  time  domain  allows  more  physi¬ 
cal  insight  into  the  processing  capability  of  the  LMS  algorithm  and  into 
its  advantages  and  limitations  for  these  kinds  of 
applications. 

Two  adaptive  deconvolution  experiments  which  used  the  storage  corre¬ 
lator  are  described.  In  the  first  experiment,  an  undesirable  time- 
delayed  ghost  pulse  is  removed.  In  the  second  experiment,  the  ringing 
in  the  pulse  response  of  an  acoustic  transducer  is  removed  by  adaptive 
filtering.  The  results  of  these  experiments  are  compared  to  computer 
simulations  with  the  LMS  algorithm.  Conclusions  are  drawn  regarding 
which  properties  and  nonlinearities  of  the  storage  correlator  are 
limiting  the  performance  of  the  adaptive  filter. 
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II.  EXPERIMENTAL  PROCEDURE 
A.  The  LMS  Algorithm 

Consider  an  adaptive  filter  with  an  input  x(t)  and  output 
y(t)  =  x  *  w  .  The  tap  weights  w  are  adjusted  after  each  iteration 
such  that  the  output  y(t)  converges  to  a  desired  signal  d(t)  .  The 
time  it  takes  for  convergence  to  occur  is  commonly  called  the  learning 
or  training  time.  After  the  optimum  set  of  weights  w  is  determined, 
the  filter  can  be  used  to  remove  the  distortions  in  signals.  For 
example,  in  a  TV  system  with  a  ghost,  the  filter  can  be  trained  on  the 
sync  signal  (which  occurs  at  the  end  of  each  line)  to  remove  the  echo 
sync  signal,  and  then  the  entire  TV  line  can  be  passed  through  the 
filter,  and  the  "ghosts"  will  be  removed. 

The  error  e(t)  between  the  desired  signal  d(t)  and  the  output  of 
the  filter  y(t)  is 


e(t)  =  d(t)  -  y(t) 


(1) 


where 


y  =  x  *  w 


(2) 


and  *  indicates  convolution.  Then  the  LMS  algorithm  specifies  that 
to  minimize  the  error,  w  must  be  changed  by  Aw  such  that 


Aw  =  )ix-A  e 


\3) 


where  indicates  correlation  and  p  is  a  constant  which  determines 

the  rate  of  convergence  of  the  system.  j 

:  1 
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Thus  it  follows  from  Eqs.  (1),  (2),  and  (3)  that 


Aw  =  yx  ■&' (d  -  x  *  w)  (4) 

A  proof  is  given  in  the  following  section  that  if  these  adjustments 
Aw  are  made,  then  the  weights  exponentially  approach  the  Wiener  filter 
solution;  this  is  a  filter  which  minimizes  the  average  mean  square  error, 
over  the  frequency  band  of  interest,  between  its  output  signal  and  a 
desired  signal  D(w)  . 

It  can  be  seen  from  the  expression  for  the  tap  weight  adjustment 
(Eq.  (4))  that  for  each  iteration,  a  convolution  with  x  ,  and  a  corre¬ 
lation  with  x  must  be  performed.  In  addition,  an  analog  storage  regis¬ 
ter  is  needed  for  w  ,  and  it  must  be  possible  to  adjust  each  of  the  tap 
weights. 

The  ASW  storage  correlator  is  ideally  suited  for  this  purpose,  as  it 
can  perform  all  of  these  operations  (convolution,  correlation,  addition, 
and  storage).  A  big  advantage  of  this  approach  is  its  simplicity.  All 
of  the  necessary  operations  are  performed  inside  a  two  (or  three)  port 
device  in  real  time. 

B.  The  Storage  Correlator 

The  ASW  storage  correlator  is  a  device  which  utlizes  the  interaction 
between  a  propagating  acoustic  surface  wave  and  a  charge  distribution 
stored  in  an  array  of  diodes  to  store,  correlate,  or  convolve  two  broad 
bandwidth  signals  in  real  time.  A  series  of  signals  may  be  added 
together  by  successively  storing  them  in  the  diode  array.  A  brief  descrip 
tion  of  the  storage  correlator  follows.  Complete  descriptions  of  the 
device*^’  **  and  its  operation  as  well  as  a  theory  of  the  device*"* 
are  given  in  the  literature. 
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A  schematic  drawing  of  a  storate  correlator  is  shown  in  Fig.  1.  If  a 

modulated  carrier  x  is  applied  to  the  acoustic  port  (#2)  and  a  signal  e 
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to  port  #3  ,  then  the  charge  distribution 

Aw  =  x  ^  e  (5) 

is  added  to  the  charge  distribution  w  stored  in  the  diode  array.  If 

a  signal  x  is  then  applied  to  port  #3,  then  the  convolution 

y  =  x  *  w  (6) 

is  the  output  at  port  #2.  The  LMS  algorithm  can  be  realized  by  repeating 
the  process  many  times.  Eventually,  the  tap  weights  w  converge  to  an 
approximation  to  the  Wiener  solution,  and  the  programming  is  complete. 

The  filter  may  then  be  used  to  deconvolve  other  signals  which  have  been 
distorted  in  the  same  manner. 

The  external  connections  used  to  operate  the  storage  correlator 
are  indicated  in  Fig.  2a.  A  much  simpler  method  is  indicated  in  Fig.  2b 
when  two  acoustic  ports  of  the  correlator  are  used,  and  the  external 
delay  line  and  the  switches  are  not  needed.  A  difference  amplifier. is  not 
needed  to  subtract  d(t)  and  y(t)  ;  rather  the  carrier  frequency  is 
adjusted  so  that  there  was  a  phase  difference  of  ir  between  the  delayed 
signal  y(t)  and  the  desired  signal  d(t)  .  This  is  explained  more 
fully  in  the  next  section. 

The  storage  correlator  is  not  an  exact  implementation  of  the  Widrow 
type  LMS  adaptive  filter  for  two  reasons.  First,  the  LMS  algorithm 
requires  that  each  tap  be  changed  individually  and  the  error  signal  for 
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the  next  tap  be  calculated  after  the  previous  tap  is  changed.  .  However, 
in  our  implementation,  the  entire  error  signal  is  calculated,  and  the 
entire  set  of  taps  is  changed  at  once.  This  difference  is  not  signifi¬ 
cant  except  for  the  case  of  very  high  gain  and  very  fast  convergence. 

The  second  difference  is  that  the  effect  of  the  plate  signal  (port  #3) 
on  the  tap  weights  is  not  linear  (Fig.  3).  If  the  error  signal  is  less 
than  a  threshold  determined  by  the  characteristics  of  the  devices  (4  V), 
the  tap  weights  are  unaffected.  The  error  signal  is  clipped  around  9V, 
so  only  ~6  dB  of  dynamic  range  is  available  at  port  #3.  If  the  error 
signal  is  fed  into  an  acoustic  port  (as  in  Fig.  4b),  then  at  least  35  dB 
dynamic  range  should  be  available. 

The  weight  signal  w(z)  is  stored  in  2000  diodes;  however,  the 
number  of  equivalent  taps  is  much  less  since  the  bandwidth  is  8  MHz 
and  the  maximum  signal  duration  is  3  ys  .  Consequently,  the  time  band¬ 
width  product  is  24.  The  device  was  represented  in  the  computer 
simulation  by  a  24  tap  transversal  filter.  The  necessary  correlations 
and  convolutions  were  calculated  in  the  computer,  and  the  error  signal 
was  modified  according  to  the  transfer  characteristic  shown  in  Fig.  3. 

It  can  be  seen  that  this  transfer  characteristic  closely  approximates 
the  experimentally  observed  characteristic. 

III.  THEORY 
A.  Introduction 

14  15 

The  LMS  algorithm  has  been  analyzed  in  detail  by  Widrow  *  .  The 

first  task  in  his  analysis  is  to  transform  from  the  coordinate  system 
consisting  of  N  samples  from  the  waveform  to  a  set  of  normal  coordi¬ 
nates  where  the  correlation  matrix  is  diagonal;  the  eigenvectors  of 
the  correlation  matrix  are  mutually  orthogonal,  and,  most  importantly. 
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the  decay  modes  are  uncoupled.  Once  this  transformation  is  accomplished, 
expressions  for  the  decay  constants  ,  misadjustment,  etc.  may  be 
derived . 

In  the  analysis  presented  here,  the  transformation  is  very  easily 
accomplished  by  taking  the  Fourier  transform  of  all  quantities.  Corre¬ 
lations  are  then  simple  products,  and  the  convergence  of  a  quantity  at 
a  given  frequency  is  independent  of  other  frequencies,  i.e.,  the  modes 
are  uncoupled. 


B.  Convergence 

Consider  the  situation  in  Fig.  2,  where  all  signals  are  rf  modulated, 
and  the  carrier  frequency  is  Wq  .  The  desired  signal  d  (t)  and  the  dis¬ 
torted  signal  xn(t)  are  repetitive,  and  as  the  number  of  iterations  n 
increases,  the  only  change  in  these  signals  is  a  phase  change 

jco  nT 

dn(t)  =  d(t)  e  U  (7) 


where  T  is  the  iteration  time.  In  the  frequency  domain 


jo>0nT 

Dn(w)  =  D(ca)  e  u  (8) 


Similar  expressions  hold  for  the  distorted  signal  X(o)) 
The  filter  output  is 


Yn(w)  =  X(oOWn(a>)  e 


jco0(nT  ♦  Tfl) 


(9) 


and  the  error  signal  is 


Equation  (14)  converges  only  if  } a{  <  1  .  The  fastest  convergence 
is  obtained  for  the  carrier  frequencies  where  the  error  signal  is  the 
subtraction  of  D  and  Y  ,  i.e.  when 


7t 


T  - 
lD 


(2p  +  1) 


(16) 


where  p  is  a  large  integer.  If  the  feedback  gain  y  is  sufficiently 
small 


1 

y  <  —  (17) 

XX* 

hen  the  convergence  requirement  }  A.{  <  1  is  satisfied  for  carrier 
frequencies  in  bands  of  width  Ao>  around  the  frequencies  given  by 
Equation  (16),  where 

2 

Aco  =  - arccos  (yXX*)  (18) 

Itd  -  To' 

It  can  be  seen  from  Equation  (17)  that  when  y  is  near  its  upper  limit, 
the  carrier  frequency  must  be  close  to  one  of  the  values  given  by 
Equation  (16). 

If  the  feedback  gain  is  significantly  less  than  its  maximum  value,  • 


it  follows  that 


JVT0  ” 

=  1  +  2yXX*  e  u  u  u 


»  e 


■1/t 


(19) 


where 


1/T  =  2pXX*  e 


jI“o(To  "  V  + 


(20) 


The  com'ergence  of  the  weight  vector  (Eq.  (14))  is  approximately  exponen¬ 
tial.  So  we  can  write 


W  =  (W,  -  W  )  e'n/T  +  W 

n  "  ^  1  w'  w 


(21) 


The  error  is 


E_  *  D  -  X*W 
n  n 


En  -  x*<Wi  "  V  e 


-n/t 


(22) 


The  error  decays  exponentially  to  zero.  In  practice,  the  bandwidth  of 
the  system  is  limited  by  the  bandwidth  B  of  the  acoustic  transducers. 
Consequently,  the  error  outside  this  frequency  range  is  not  affected  by 
the  adaptation  process.  It  will  also  be  noted  that  phase  distortion  in 
-the  system  is  cancelled  out  by  the  application  of  the  convolution  and 
correlation  process  in  turn.  So,  although  this  is  basically  a  feedback 
system,  there  is  no  problem  with  instability  provided  the  convergence 
criteria  are  satisfied. 

We  have  implicitly  assumed  that  the  signals  x(t)  and  d(t)  are 
sufficiently  short  compared  to  the  time  length  of  the  diode  array  that 
all  of  the  signal  M>(t)  is  stored  in  the  diode  array.  If  this  is  not 

-  11  - 


the  case,  then  different  frequencies  do  not  decay  independently,  and 
the  analysis  is  more  difficult.  An  important  result  of  this  more 
difficult  case  is  that  the  error  does  not  decay  to  zero  in  the  passband. 

If  noise  on  the  input  signal  X  is  included  in  this  analysis,  and 
if  the  noise  on  the  plate  and  acoustic  ports  are  correlated,  as  they 
would  be  if  a  delay  line  were  used  to  generate  the  second  X  signal 
which  is  needed  during  each  iteration,  then  the  expression  for  the  con¬ 
vergence  of  the  weight  signal  is  again  given  by  Equation  (22)  where 

DX* 

w  =  -  (23) 

w  XX*  +  <NN*> 

1 

T  =  -  (24) 

2p  (XX*  +  <NN*>) 

Equation  (23)  is,  of  course,  the  Wiener  filter  solution. 

IV.  RESULTS  AND  DISCUSSION 
A.  Echo  Suppression 

In  this  experiment,  a  square  pulse  0.4  ps  long  is  followed  by  an 
echo  pulse.  The  desired  signal  is  a  single  pulse  0.4  ps  long.  For 
an  echo  which  is  6  dB  less  in  amplitude  than  the  main  pulse  (Fig.  4a), 
the  sidelobe  suppression  after  10  iterations  (200  ps)  is  15  dB,  as 
shown  in  Fig.  4b.  The  dependence  of  sidelobe  suppression  on  echo  height 
is  shown  in  Fig.  5.  The  results  of  computer  simulations  of  the  LMS 
algorithm  with  22  taps  and  using  clipping  and  a  threshold  6  dB  below 
the  clipping  level  are  also  shorn  in  Fig.  5.  The  computer  simulation 
agrees  very  well  with  experimental  results  except  that  the  maximum  side¬ 
lobe  suppression  is  23.5  dB  which  is  4.5  dB  higher  than  was  experimentally 
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obtained.  The  reason  is  that  the  spurious  signals  generated  in  the  device 
during  read  out  limits  the  dynamic  range  of  the  signal  input. 


An  important  result  obtained  in  this  experiment  is  that  spurious 
acoustic  signals  generated  by  the  plate  readout  signal  can  be  suppressed 
by  up  to  13  dB  as  a  result  of  the  adapting  process.  The  filter  does 
not  distinguish  between  echoes  and  distortions  generated  externally  or 
by  the  device  itself.  This  result  is  demonstrated  in  Fig.  6.  The  upper 
trace  is  the  adaptive  filter  result  after  removing  an  echo  from  a  0.4  ys 
long  pulse.  If  all  signals  except  the  plate  readout  signal  are  removed, 
the  filter  output  (lower  trace)  is  the  spurious  signal  generated  by  the 
plate  signal.  A  large  spurious  signal  can  be  seen  when  previously  there 
had  been  a  null. 

The  advantage  of  computer  simulation  of  this  adaptive  filter  is  that 
the  threshold  and  clipping  levels  may  be  easily  changed  to  see  what  effect 
they  have  on  the  performance  of  the  adaptive  filter.  These  results  are 
summarized  in  Table  I.  The  computer  simulated  linear  LMS  result  is 
given  in  the  first  row.  If  a  threshold  level  is  included,  then  the 
computer  simulation  converges  much  faster,  but  to  worse  results.  The 
values  obtained  are  in  agreement  with  experimental  results.  If  the 
feedback  gain  is  increased  so  that  the  error  signal  is  now  clipped, 
then  much  better  sidelobe  suppression  is  obtained  in  both  the  experimental 
and  computer  simulation  cases.  We  note  that  clipping  of  the  signal 
increases  the  rate  of  convergence  radically,  as  has  been  noted  by  others. 
The  algorithm  employed  is  therefore  known  as  the  clipped  LMS  algorithm. 

Regardless  of  the  shape  of  the  desired  signal  (d)  and  input  signal 
(x),  it  was  always  experimentally  observed  that  the  feedback  gain  must  be 
large  enough  to  strongly  clip  the  error  signal  during  the  first  few 
iterations  for  the  optimum  filtering.  After  many  iterations,  the  error 
signal  is  only  slightly  clipped. 
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The  threshold  and  clipping  levels  can  be  changed  by  modifying  the 
design  of  the  storage  correlator.  For  this  reason,  a  series  of  computer 
simulations  was  made  to  determine  the  effect  of  having  other  threshold 
to  clipping  ratios  (o)  than  the  value  a  -  .5  which  exists  in  the 
present  devices.  Only  the  ratio  a  is  important  here.  The  minimum 
r.m.s.  error  does  not  change  if  the  threshold,  clipping,  and  gain  values 
are  all  multiplied  by  the  same  value.  The  asymptotic  value  of  the  r.m.s. 
error  is  shown  as  a  function  of  a  in  Fig.  7  for  the  case  of  a  6  dB 
echo.  The  optimum  value  for  the  gain  is  used  at  each  point.  The  end 
point  o  =  1.0  corresponds  to  the  situation  of  hard  clipping  where  the 
threshold  level  equals  the  clipping  level,  and  the  feedback  error  is 
either  1  or  0.  The  r.m.s.  error  increases  linearly  with  a  .  This  is 
intuitively  expected  since  any  error  values  less  than  the  threshold 
do  not  affect  the  tap  weights  and  are  not  adapted  to  zero.  The  r.m.s. 
error  would  then  linearly  increase  with  threshold  level.  Consequently,  the 
design  of  the  storage  correlator  should  minimize  the  threshold  value. 

B.  Reduction  of  Bulk  Transducer  Ringing 

The  object  of  this  experiment  was  to  improve  the  impulse  response  of 
an  acoustic  bulk  wave  1.25  MHz  transducer  used  for  acoustic  nondestructive 
evaluation.  The  impulse  response  is  shown  in  Fig.  8a.  The  desired  signal 
is  a  unipolar  pulse  with  a  width  of  .3  ys  .  If  the  impulse  response  is 
correlated  with  itself  (Fig.  8b),  the  highest  sidelobe  is  1.7  dB  smaller 
than  the  peak.  However,  if  the  correlator  is  used  as  an  adaptive  filter, 
then  the  sidelobe  level  can  be  reduced  to  7  dB  below  the  peak  after  10 
iterations  (Fig.  8c)  and  10  dB  below  the  peak  after  35  iterations. 

The  growth  of  the  peak  as  the  number  of  iterations  increases  is 
shown  in  Fig.  9.  The  computer  simulation  is  also  shown  and  can  be  seen 
to  be  in  good  agreement  with  the  experimental  result.  The  computer 


simulation  is  in  good  agreement  with  the  experimental  result  except 
that  the  maximum  predicted  suppression  is  2.6  dB  better  than  was  experi¬ 
mentally  observed. 

Linear  LMS  theory  predicts*  that  the  time  averaged  mean  square 
error  (MSE)  should  decay  exponentially  to  a  constant  value.  The  ratio 
of  this  constant  value  to  the  Wiener  solution  is  called  the  misadjustment 
M  .  The  misadjustment  is  approximately  related  to  the  decay  constant  T 
by  the  relation* 


M  =  N/4t  (2s) 

where  N  is  the  number  of  tap  weights.  In  our  case,  N  is  the  time 
bandwidth  product  of  the  storage  correlator.  Thus,  the  faster  the 
convergence,  the  larger  the  misadjustment. 

The  experimentally  observed  decay  of  the  MSE  for  the  bulk  transducer 
case  is  shown  in  Fig.  9.  The  decay  is  exponential  except  for  the  first 
few  iterations  when  the  error  signal  is  strongly  clipped  and  the  mean 
square  error  decays  faster  than  an  exponential.  The  experimental  decay 
constant  is  7.5,  and  the  misadjustment  calculated  from  Equation  (25)  is 
equal  to  .S  .  We  interpret  this  to  mean  that  the  final  result  is  close 
to  the  Wiener  solution. 

C.  Limitations 

The  bandwidth  of  this  device  (8  MHz)  is  larger  than  the  bandwidth  of 
other  methods  of  implementing  the  LMS  algorithm;  however,  larger  bandwidths 
are  often  desirable.  Our  technique  can  work  well  as  an  adaptive  filter 
for  a  1.25  MHz  transducer,  but  it  works  poorly  for  a  2.5  MHz  transducer. 

The  diode  array  length  (3  ps)  is  also  a  limitation  since  we  cannot 
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adapt  signals  which  are  longer  than  this.  Also,  since  the  correlation 
signal  is  truncated,  the  filter  output  is  distorted.  These  two  limi¬ 
tations  are  manifestations  of  the  fact  that  if  the  correlator  had  a 
larger  time  bandwidth  product  than  24,  then  it  would  equivalently  have 
more  taps  and  could  adapt  a  broader  class  of  signals.  Time  bandwidth 
products  of  500  or  more  are  possible  with  the  use  of  present  technology. 
Note,  however,  that  Equation  (25)  predicts  that  if  the  number  of  taps 
is  increased,  then  the  convergence  time  also  increases  (for  a  given 
level  of  misadjustment) .  Thus,  a  larger  number  of  taps  is  not  desirable 
for  applications  in  which  the  distortion  is  changing  rapidly.  For 
example,  if  the  timing  or  amplitude  of  the  echo  significantly  changes 
over  a  time  period  of  200  ys  (10  iterations),  then  a  larger  number  of 
taps  would  not  be  desirable. 

V.  CONCLUSIONS 

It  was  demonstrated  that  a  -6  dB  echo  could  be  suppresed  to  -15  dB 
in  200  ys  (10  iterations).  The  sidelobe  of  the  impulse  response  of  a 
1.25  MHz  acoustic  bulk  transducer  can  be  reduced  from  -2  dB  to  -10  dB 
with  adaptive  filtering. 

The  experimental  results  are  on  the  whole  in  excellent  agreement  with 
the  computer  simulations.  Thus  it  appears  that  the  storage  correlator 
is  being  operated  in  the  optimum  manner  and  that  the  only  major  distor¬ 
tion  introduced  by  the  storage  correlator  is  the  presence  of  the  threshold 
and  clipping  levels. 

There  are  many  advantages  of  using  a  storage  correlator  to  implement 
the  clipped  LMS  algorithm  for  adaptive  filtering.  Only  two  connections 
are  needed  for  a  device  which  can  correlate,  convolve,  store,  and  add 


2 

broadband  rf  signals  in  real  time.  The  device  is  small  (1  cm  )  and 
has  the  potential  of  much  broader  bandwidths  (100  MHz)  than  is  obtain¬ 
able  with  digital  techniques.  The  device  is  far  faster  in  operation 
than  any  competitive  inverse  filter  with  its  bandwidth. 
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TABLE  I 


Experimental  Results  and  Computer  Simulations  for  a 
Filter  Input  of  a  .4  ys  Wide  Pulse  with  a  -6  dB 
Echo.  The  Computer  Simulation  Uses  24  Taps. 
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FIGURE  LEGENDS 


1.  Schematic  drawing  of  storage  correlator. 

2.  (a)  Actual  and  (b)  proposed  implementation  of  SAW  adaptive  filter. 

3.  Correlator  output  as  a  function  of  plate  (readout)  pulse  height. 

4.  (a)  Input  signal  and  (b)  output  signal  after  22  iterations. 

5.  Experimental  and  computer  simulation  results  for  echo  suppression 
experiment . 

6.  Echo  reduction  after  10  iterations  (upper  trace).  Spurious  signals 
generated  by  plate  read-out  signal  (lower  trace). 

7.  Minimum  r.m.s.  error  as  a  function  of  threshold  level. 

8.  (a)  Impulse  response  of  1.25  MHz  bulk  acoustic  transducer;  (b)  auto¬ 
correlation  of  impulse  response;  and  (c)  filter  output  after  10 
iterations  of  adapting. 

9.  Output  pulse  level  during  adapting  process. 

10.  Mean  Square  Error  during  adapting  process. 
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